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Abstract: The Australian Square Kilometre Array Pathfinder (ASKAP) will give us an unprecedented 
opportunity to investigate the transient sky at radio wavelengths. In this paper we present VAST, an 
ASKAP survey for Variables and Slow Transients. VAST will exploit the wide-field survey capabilities 
of ASKAP to enable the discovery and investigation of variable and transient phenomena from the 
local to the cosmological, including flare stars, intermittent pulsars. X-ray binaries, magnetars, extreme 
scattering events, interstellar scintillation, radio supernovae and orphan afterglows of gamma ray bursts. 
In addition, it will allow us to probe unexplored regions of parameter space where new classes of 
transient sources may be detected. In this paper we review the known radio transient and variable 
populations and the current results from blind radio surveys. We outline a comprehensive program 
based on a multi-tiered survey strategy to characterise the radio transient sky through detection and 
monitoring of transient and variable sources on the ASKAP imaging timescales of five seconds and 
greater. We also present an analysis of the expected source populations that we will be able to detect 
with VAST. 

Keywords: telescopes — surveys — radio continuum: general — galaxies: general — stars: general 
— ISM: general 
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1 Introduction 

Many classes of astronomical objects are known to be 
variable radio sources, including the Sun, the planets, 
cool stars, stellar binary systems, neutron stars, su- 
pernovae, gamma-ray bursts and active galactic nuclei. 
There have been extensive studies of these objects us- 
ing targeted surveys, but few blind, unbiased surveys 
for variable radio phenomena. An effective survey for 
radio variability needs to maximise the metric 



Table 1: ASKAP Specifications 



(1) 



where A is the effective collecting area of the telescope, 
n is the solid angle coverage of the survey, T is the total 
dura tion of the observations and At i s the time resolu- 
tion (|Cordes et al.ll2004l : [007569120071 ). Typically blind 
radio surveys have had either too small a field of view 
and too little time coverage per field of view to al- 
low a comprehensive survey for transient and variable 
phenomena. 

The Square Kilome tre Array (SKA n is a proposed 
future radio telescope (jPewdnev et al.lEoOgI ') that wiU 
be many times more sensitive than any existing facil- 
ity. Should it go ahead, the SKA will be designed to 
explore five key science projects: galaxy evolution, cos- 
mology and dark energy; strong-field tests of gravity 
using pulsars and black holes; the origin and evolu- 
tion of cosmic magnetism; probing the dark ages — 
the first black holes and stars; and the cradle of life, 
searching for life and planets. In addition, there is 
an awareness of building in (or not designing out) the 
possibility of serendipitous discoveries. In particular, 
the dynamic radio sky is recognise d as a rich and rela- 
tively unexpl ored discovery space (|Cordes et al. I l2004l : 
ICorded[2007h . Plans for the SKA are still under de- 
velopment, with a Phase 1 instrument expected to be 
completed in 2020, and a Phase 2 instrument in 2024. 

The Australian Square Kilometre Array Pathfinder 
(ASKAP J3 

is a precursor and technology development 
platform for the full SKA and will greatly improve 
on our ability to detect radio transients and variable 
sources. Its moderately high sensitivity (RMS of ~1 mjy 
in 10 seconds), combined with a wide field of view (cov- 
ering 30 square degrees of sky in a sing le pointing) will 
enable fast, sensitive all-sky surveys (|johnston et al.l 
[2007 . 2008). ASKAP is currently under construction 
in Western Australia. The final telescope will have 
36 12-metre antennas, with a maximum baseline of 
6 km. The wide field of view is enabled by the use 
of Phased Array Feeds (|Chippendale et aITl2010l ) of 
100 dual-polarisation pixels that will operate over the 
range 700-1800 MHz. A summary of the planned 
ASKAP specifications is given in Table [T] The 6- 
antenna Boolardy Engineering Test Array (BETA), is 
expected to be completed in late 2012 and will allow 
some early testing of the techniques required for the 
full ASKAP surveys. 

ASKAP is primarily a survey instrument, and this 
is refiected in the time allocation policy. In its initial 



Number of dishes 


36 


Dish diameter (m) 


12 


Dish area (sq m) 


113 


Total collecting area (sq m) 


4072 


Aperture efficiency 


0.8 


System temperature (K) 


50 


Field-of-view (deg^) 


30 


Frequency range (MHz) 


700-1800 


Bandwidth (MHz) 


300 


Maximum number of channels 


16 384 


Maximum baseline (km) 


6 


RMS cont. sensitivity (10 s) 


640/.iJy/bm 


RMS cont. sensitivity (1 hr) 


47^Jy/bm 



five years of operation 75% of the available time will be 
scheduled for several major Survey Science Project^ 
that w ere selected through a competitive call for pro- 
posals (|Ballll2009l ). 

In this paper we describe one of these projects. 
An ASKAP Survey for Variables and Slow Transients 
(VASTfl In the context of ASKAP, 'slow' means vari- 
ability or transient behaviour detectable on timescales 
greater than the cadence on which images will be per- 
formed (expected to be 5-10 seconds) . A related project. 
The Commensal Real-Time ASKAP Fa st-Transients 
Survey (CRAFT; iMacguart et al.ll2010l ). will search 
for 'fast' transients on timescales shorter than 5 sec- 
onds. Different techniques are required to process and 
analyse data for fast and slow transients, and from a 
scientific perspective this division roughly represents a 
separation between coherent and incoherent emission 
processes. The main science goals of VAST are: 

• to detect and monitor 'orphan' gamma-ray burst 
afterglows to understand their nature; 

• to conduct an unbiased survey of radio super- 
novae in the local Universe; 

• to determine the origin and nature of the struc- 
beam tures responsible for extreme scattering events; 

• to make a direct detection of baryons in the in- 
tergalactic medium; 

• to discover flaring magnetars, intermittent or 
deeply nulling radio pulsars, and rotating ra- 
dio transients through changes in their pulse- 
averaged emission; 

• to investigate intrinsic variability of active galac- 
tic nuclei; 

• to detect and monitor flare stars, cataclysmic 
variables and X-ray binaries in our Galaxy; and 

• to discover previously unknown classes of ob- 
jects. 

To achieve these goals VAST will conduct three sur- 
veys, in addition to commensal observing with other 



-"^See http://www.skatclcscope.orgj 

^See jhttp: / /www.atnf.csiro.au/projects/askap[ 



•^Thc full list of projects is availab le at 
|http: / / www.atnf.csiro.au/projects/a skap/ssps.htmll 
''See |http://www. askap.org/vast I 
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Table 2: Survey parameters for the VAST surveys (see Section HTTI for details). The two components of 
the VAST-Decp survey are discussed further in the text. The Commensal column gives an indication 
of the amount of time that VAST will have access to through commensal observing with other ASKAP 
projects. 





VAST-Wide 


VAST-Deep 
Multi-field Single field 


VAST-Galactic 


Commensal 


Observing time (hrs) 


4380 


3200 


400 


600 


1.5 years 


Survey area (sq deg) 


10 000 


10 000 


30 


750 


10 000 


Time per field 


40 s 


1 hr 


1 hr 


16 min 


12 hours 


Repeat 


daily 


7 times 


daily 


64 times 


none 


Observing freq (MHz) 






1130-1430 






Bandwidth (MHz) 






300 






RMS sensitivity (mJy beam~^) 


0.5 


0.05 




0.1 


0.01 


Field of view (sq deg) 






30 






Angular resolution 






10" 






Spectral resolution 






10 MHz 






Time resolution 






5 seconds 






Polarisation products 






IQUV 







ASKAP projects: VAST-Wide, aimed at detecting rare, 
bright events, will survey 10 000 square degrees per 
day to an RMS sensitivity of 0.5 mJy beam~^; VAST- 
Deep, aimed at detecting unknown source classes, will 
survey 10 000 square degrees down to an RMS sensitiv- 
ity of 50/iJy beam"^; and VAST-Galactic which will 
cover 750 square degrees of the Galactic plane, plus 
the Large and Small Magellanic Clouds down to an 
RMS sensitivity of 0.1 mJy beam"^. The full survey 
parameters are shown in Tableland discussed in Sec- 
tion g] 

In Section[2]we motivate the VAST surveys by dis- 
cussing the science goals in more detail, and explain 
how they fit in the context of existing efforts. We also 
make predictions about the likely source populations 
that will be observed by VAST. In Section[3]we review 
previous and current blind radio surveys, as well as 
archival studies. In Section|4]we describe the proposed 
VAST survey strategy in detail. Finally, in Section [S] 
we discuss the design of the VAST transient detection 
pipeline. 



2 VAST Science 

Astronomical transient phenomena are diverse in na- 
ture, but they can be broadly classified on the ba- 
sis of their underlying physical mechanism into one 
of four general categories: explosions, propagation ef- 
fects, accretion-driven, and magnetic field-driven events. 
In this section we discuss the main VAST science goals 
in each of these areas and discuss the expected popu- 
lations for different types of sources. 



2.1 Explosions 

Gamma-ray bursts (GRBs) and supernovae (SNe) are 
some of the most energetic explosions in the Universe. 
Although they have been studied extensively there are 



a number of unresolved questions about their na ture, 



inclu ding the existence of radio orphan afterglows (iLevinson et al.l 
2002 ) , the beaming frac tion of gamma-ray bursts (|Sari et al.l 
199S : iDalal et al.l 20021 ) . and the current rate of mas- 
sive star formation in the Universe as traced by the 
cosmological radio supernova rate. These questions 
can be addressed by surveys that observe a large area 
of sky to good sensitivity, on a regular basis. The ca- 
pacity of VAST to detect these rare objects is orders of 
magnitu de greater than previo us blind surveys such a s 
those of I Gal- Yam et all (|2006l 'l or lBower etal] (|2007f ). 



2.1.1 Gamma- Ray Bursts 

The afterglow emission from gamma-ray bursts is gen- 
erated via particle acceleration in the decelerating blast 
wave. At late times, the peak of the afterglow spec- 
trum is at radio energies, and thus radio studies pro- 
vide the most interesting probes for GRBs scie nce. 
Based on previous studies (|Chandra fc Frailll2012l ). we 
expect that the GRBs seen by ASKAP will be long- 
lasting but intrinsically fainter on average than shorter 
wavelength afterglows. Given the sensitivity of ASKAP, 
this means that those GRBs studied by VAST will, in 
general, be the rare but important nearby (~1500 Mpc) 
events. In the relativistic expansion phase (up to 100 
days after the initial explosion), ASKAP monitoring in 
conjunction with optical and X-ray facilities will allow 
us to measure physical parameters such as the total 
kinetic energy of the explosion, the density of the cir- 
cumburst medium and the geometry of the outflow. 

At late times, when the blast wave has deceler- 
ated to sub-relativistic velocities, the jetted outflows 
expand and become more isotropic. Under these cir- 
cumstances, the the total (i.e. calorimetric) energy of 
GRB blast waves may be e stimated - a key but diffi- 
cult parameter to measure (|Frail et al.ll2005l ). Instead 
of relying on rare, nearby events, VAST may also ob- 
tain a robust beaming-independent estimates GRB en- 
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ergies from a large sample of single-epoch flux density 
measurements (and limit s) of afterglows at late times 
(|Shivvers fc BergeJ[2oTll ). 

As the relativistic beaming is totally suppressed 
at late time, radio monitoring is the best solution to 
identify 'orphan afterglows', i. e., those GRB s whose 
jet is beamed away from us (|Rhoadsl Il997l '). While 
they would lack a gamma-ray trigger by definition, 
such systems could be distinguished from other classes 
of transients through the evolution and luminosity of 
their light curves, their location towards galaxies that 
resemble other GRB hos ts, and their probable asso- 
ciation with supernovae (|Soderberg et al.ll2010l ). The 
discovery of a population of orphan afterglows would 
allow us to constrain the beaming fraction of GRBs - 
a key quantity for understanding the true GRB event 
rate. It would also open up an order of magnitude 
more objects which can be scrutinised in detail. This 
may also give insights on the jet structure of GRBs and 
clarify t he puzzling lack of j et breaks from the Swift 
satellite (|Racusin et all 120091 '). 

2.1.2 Supernovae 

While thermonuclear (Type la) supernovae have proven 
to be extremely useful in measuring the acceleration of 
the Universe, no T ype la SN has ever been detected at 
radio wavelengths (iPanagia et al.ll2006l : [Hancock et al.l 
l201ll : IChomiuk et al.ll2012l ). indicating a very low den- 
sity of circumstellar material. In contrast, many core- 
collapse SNe (those of Type lb, Ic, and the various 
Type II sub-classes) h ave been detected and monitored 
at radio wavelengths (|Weiler et al.l |2002| : iRvder et al.l 
l2004l : [SoderberdBoOTl ). Radio observations replay the 
late phases of stellar evolution (in reverse), and reveal 
clues about the late-time evolution of the progenitor 
system. 

With VAST we will conduct an unbiased census 
of core-collapse SNe (a still undetermined fraction of 
which emit at radio wavelengths), allowing us to match 
radio detections against optically-discovered SNe. We 
will also detect new radio SNe that may have gone un- 
detected in the o ptical or infrared due to significant 
amounts of dust (|Kankare et al.l |2008| ) . The preva- 
lence of such heavily-obscured SNe is virtually uncon- 
strained, with hints that they may constit ute a non- 
trivia l subpopulation. Recently, SN 2008cs (|Kankare et al.l 
l2008l ) was detected in the IR and con firmed in the 
radio bands, and SN 2008iz (Brunthal cr et all |2009| ) 
was serendipitously discovered in radio observations 
of M82. SN 2008iz has remained undetected in the 
optical, IR, and X-ray b ands despite sensitive searches 
(jBrunthaler et al]l2010l ). Discovery of new SNe is valu- 
able for tying down the current rate of massive star 
formation in the Universe and for understanding how 
and when SNe explode in star-bursting galaxies. Radio 
emission from core-collapse SNe is detectable within 
days, and they reach their peak luminosity at 1.4 GHz 
as much as one year after the explosion, making them 
ideal targets to detect and study with the sensitivity 
and cadence of VAST (see Section |4]). 

Beyond the large population of standard core-collapse 
SNe, there exist sub-classes of objects that are of par- 



ticular interest. For instance, radio observations of a 
recent type Ib/c supernova, SN2009bb, showed that 
it was expanding at mildly relativisti c velocities (/3 ~ 
0.8) powered by the central engine (jSoderberg et al.) 
I2OIOI ). Such SNe are usually associated with GRBs and 
are discovered by their intense but short-lived gamma- 
ray energy emission. Yet, SN2009bb had no detected 
gamma-ray counterpart. It has been suggested that 
the blast-waves produced in such engine-driven SNe 
could be responsible for the acceleration of u ltra-high 
energy cosmic-rays (jChakraborti et aP I2OIII ) . How- 
ever, the present limits on the neutrino flux expected 
as a by-product of cosmic-ray productio n in GRBs are 
at lea st a factor of 3.7 below predictions (|Abbasi et al.l 
l2012i ). Radio surveys with the capability of complet- 
ing the sample of these explosions will improve our un- 
derstanding of particle acceleration processes and the 
diversity of core-collapse events (see lLien et al.|[201ll , 
for predictions on the detection rates of core-collapse 
radio SNe for a range of SKA pathfinder instruments) . 

2.1.3 Gravitational Wave Counterparts 

Gravitational wave (GW) events are expected to be 
compact and involve significant energies — both GRBs 
and SNe are considered to be among likely sources of 
both electromagnetic and GW emission. In an analo- 
gous manner, other classes of GW events may also gen- 
erate significant transient radio emission, especially if 
they occur in dense environments. We discuss the par- 
ticular case of mergers of neutron star binaries in Sec- 
tion l2.3.?1 For all classes of GW events, however, local- 
isation will be a significant challenge, with search areas 
expected t o be of order 10 -100 deg^ for planned GW 
detectors (|Fairhurstl I2OIII ). The large instantaneous 
sky coverage of ASKAP makes it particularly valuable 
in the search for electromagnetic counterparts to GW 
events, although there are significant uncertainties in 
the requisite search sensitivities and timescales. The 
search area estimates me ntioned above are computed 
for the Advance d LIGO (iHarrv et al.|[2oTol ) and Ad- 
vanced VIRGO (|Accadia et al.ll201ll ) detectors. These 
systems are planned to come online in the 2014-2015 
timeframe, meaning that ASKAP will be well-positioned 
to make discoveries in this parameter space. 

2.2 Propagation 

Plasma Propagation Effects: At centimetre wave- 
lengths any object with an angular size < 1 mas is 
subject to the effects of interstellar scintillation (ISS), 
caused by inhomogeneities in the ionised interstellar 
medium (ISM) of our Galaxy. Refractive ISS (RISS) 
is manifested as slow intensity variations with a wide 
correlation bandwidth that is caused by large-scale re- 
fracting structures in the ISM. Much faster variations 
that have narrow frequency scales are associated with 
diffraction from small-scale irregularities (DISS). Ob- 
servations of the intensity fluctuations caused by ISS 
can be used to infer both the small-scale structure of 
the ISM and the structure of a scintillating source on 
scales down to ~ 10 — 200 microarcseconds for com- 
pact extragalactic sources and ~ 10 nanoarcseconds 
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for some pulsars. 

The density inhomogeneities in th e ionised ISM are 
typic ally characterised as a power law (|Armstrong et aP 
1 19951 ). However, a small number of pulsars have dis- 
played strong fringing events in their dynamic spec- 
tra, in which the normal, random appearance of the 
dynamic spectrum caused by DISS is modified by the 
appearance of quasi-periodic fringes. Both pulsars and 
AGN have displayed extreme scattering events (ESEs) 
in their light curves, in which their flux densities ex- 
hibit marked decreases (~ 50%) bracketed by modest 
increases; and a number of AGN have exhibited intra- 
day variability (IDV), marked by flux density varia- 
tions on time scales of hours. Fringing events are in- 
terpreted most naturally as the beating of two distinct 
images of the pulsar and both ESEs and IDV can be 
interpreted as due to lensing events through a plasma. 
All of these phenomena could occur if there is a tem- 
porary increase in the amount of power on scales of 
order 10^^ cm along or close to the the line of sight, 
e.g., as due to a discrete "cl oud" of electrons drifting in 
front of the pulsar or AGN (|Cordes fc Wolszczan|[l986l : 
IClegg et all] 19981 ). In the case of IDV, it may also be 
the case that transient, high brightness temperature 
structures are required to be produced in the cores of 
some AGN. 

The characteristics of the discrete structures in- 
ferred to be responsible for these phenomena are such 
(density > 100 cm~^ , temperature > 10^ K) that they 
would be highly overpre ssured with respect to th e bulk 
properties of the ISM (|Kulkarni fc HeilesI figsi ) . and 
would likely be transient structures. While there is 
some evidence to suggest that a small fraction of the 
volume of the ISM is at a high pressure (jjenkins fc Tripnl 
1200 ll ). the distribution of such highly overpressured 
structures is not known, although the characteristics 
of IDV in some AGN suggest that some structures are 
extremely close to the Sun (~ 1 — 10 pc). The exis- 
tence of such highly over-pressured structures in the 
local ISM is hard to reconcile with our c urrent under- 
standing of interste llar turbulence (e.g. I Fiedler et al.l 
Il987bl : Imikerll2007r ). Alternatively, these phenomena 
may indicate that a simple power law is not a good de- 
scription of the spectrum of density inhomogeneities, 
which would then constrain the nature of the plasma 
processes involved in the generation and maintenance 
of the density inhomogeneities. 

Because these propagation events are transient in 
nature, many of the previous observations have been 
"accidental," happening to catch a source displaying 
such an event during the course of observations de- 
signed for other reasons. There have been relatively 
few surveys designed speciflcally to try to find prop- 
agation events, with notable examples being the Green 

Bank Interferometer (GBI) monitoring program jFielder et a jent that ESEs are caused by refractive lens-like den 
19M; Lazio ct al. 2001l3) and the MASIV survey p7ovell et alji ty perturbations along the Uiie-of-sig ht to a back- 
20031 . l2008i ). Moreover, because of either limited sensi- groun d source (e.g. iFiedler et al.|[l987bl : iRomani et al.l 

Il987l ). The refracting lenses which cause ESEs must 
be Galactic in origin and probably lie within a few kpc, 
but in the 24 years since their discovery, no satisfactory 
physical model has emerged for the ESE phenomenon. 
Analysis of existing data suggests that ESEs are un- 
like any other identified component of the ISM and 



agation events and the likelihood of a source exhibit- 
ing such an event. In turn, these much stronger con- 
straints can be used to infer properties of the ionised 
ISM and the mechanisms by which such structures are 
generated. As an illustration of the kind of monitoring 
program that ASKAP will enable, a survey could be 
designed to produce daily flux density measurements 
of sources stronger than about 10 mjy with the capa- 
bility to detect 50% flux density changes at the lOcr 
level. At this sensitivity level, a survey of even only 
1000 deg^ would still be able to monitor nearly com- 
pact 50,000 sources. 

Gravitational Lensing: Aside from the plasma prop- 
agation effects discussed in more detail below, we also 
note that strong gravitational lensing of active galax- 
ies by galaxy clusters can lead to an achromatic time- 
domain signal that contains information about both 
the active galaxy and about the intervening space- 
time. Radio observations of strong lensing have been 
conducted successfully with existing facilities and those 
data ha ve been used to constrain cosmol ogical parame- 
ters (e.g lFassnacht et al.lll999l : lMvers etral.j,2003 '). These 
operate by spatially separating the lensed components 
of the active galaxy. In general ASKAP will not have 
the angular resolution to observe such multiple imag- 
ing, but the large number of sources monitored with 
regular cadence means that auto-corre lation analysis 
(as performed bv lBarnacka et al]|201lf ) can itself con- 
strain lens time delays. 

2.2.1 Pulse Broadening from Scattering 

The well known temporal broadening of pulsar pulses 
caused by scattering from electron-density variations 
will affect any source of short duration bursts, includ- 
ing those that will be targeted by VAST. Pulse broad- 
ening from sources in the Galactic center can be broader 
than 1000 sec at 1 GHz. Bursts from extragalactic 
sources at cosmological distances can similarly be broad- 
ened, especially if the line of sight intersects another 
galaxy but also from the intergalactic medium (IGM). 
Scattering in the Milky Way, in other galaxies, and 
in the IGM therefore define pulse-broadening horizons 
beyond which pulses with s pecific widths will not be 
seen (see, e.g.. ]Cordej|2012l ). Pulse broadening from 
scattering can be distinguished from source-intrinsic 
effects by identifying across the ASKAP band the strong 



frequency dependence 



of scatter broadening. 



2.2.2 Extreme Scattering Events 

Extreme scattering events (ESEs) are flux density vari- 
ations marked by decreases in flux density bracketed 
by increases (Figure [l}. There is widespread agree- 



tivity, limited field of view, or both, most surveys have 
only been able to monitor a relatively small number 
(~ 100) of relatively strong sources (> 100 mJy). 

A wide-area, high cadence, sensitive survey would 
be able to place strong constraints on both the spa- 
tial distribution of the structures responsible for prop- 
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-400 -200 200 400 
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Figure 1: An Extreme Scattering Event in 
Q0954+658 at 2.7 GHz (lower) and 8.1 GHz (upper) 
adapted from iFiedler et al.l (|l987bl ). For clarity, an 
offset of 1 Jy has been added to the top trace. The 
strong frequency dependence of ESEs and the neces- 
sity of regular sampling of the light curves over a long 
period is evident. At ASKAP frequencies (~ 1 GHz), 
the amplitude of the flux density decrease would likely 
have been even larger, and the flux density would likely 
have increased to an even higher value during the start 
of the event. 



that they are present in vast numbers in the Galaxy 
(with a density ~ lO'^pc"^). If the refracting proper- 
ties are taken at face value, they imply regions with 
pressures ~ 1000 times higher than the diffuse ionized 
ISM. One suggestion is that they ar e associated w ith 
neutral, self-gravitating gas clouds (|Walkeill2007l ). m 
which case they would make a substantial contribu- 
tion to the total baryonic mass content of the Galaxy. 
Alternatively, they may represent small-scale density 
enhancement s created by converging flows in a turbu- 
lent medium (| Jenkins fc TripnlbOOlf ). but it is unclear 
whether such structures could persist for the observed 
ESE durations. 

A wide-area survey, with daily cadence would al- 
low us to address the physical composition and origin 
of the objects responsible for ESEs in a comprehensive 



Population stat istics: One of the m ore comprehen- 
sive ESE surveys (|Fiedler et al.|[l994bl ') provided a to- 
tal coverage of approximately 600 source-years; VAST 
has the potential to achieve a coverage that is larger 
by a factor of roughly 300, searching for ESEs toward 
sources with much lower flux densities. In turn, in- 
creased coverage will constrain the spatial density and 
distribution of refracting lenses in the Galaxy, includ- 
ing whether th ey are associated wit h particular Galac- 
tic structures (|Fiedler et al.|[l994al 'l. A larger number 
of ESEs will also allow investigation into any poten- 
tial diversity in light curve shapes, and whether any 
such di versity can be un derstood in terms of simple 
models (|Clegg et al.ll 19981 ) or the extent to which more 
complex models of the density structures, background 
sources, or both are required. 



Real-time characterisation of ESEs: A straight- 
forward approach to elucidating the nature of the struc- 
tures responsible for ESEs would be to conduct ad- 
ditional observations of a source while it is undergo- 
ing an ESE. Events like tho se seen towards qua sars 
B0954-f658 and B1749+ 096 (Fiedler ct al."l994bl) or 
B1741-038 (|Fiedler et al . 1992; Clcgg__cl_aL 1996,) should 
be seen at the rate of roughly 60 per year in an all-sky 
monitoring survey with ~1 mjy RMS sensitivity, with 
~10 in progress at any given moment. Thus, in con- 
trast to the GBI monitoring program, for which only 
one ESE was identified while it was on-going, VAST 
will enable real-time follow-up. 

Some key follow-up measurements could be made 
at higher frequencies, particularly in the optical to 
X-ray bands, which should directly reveal the dense 
neutral gas clouds responsi ble for ESEs v ia scattering, 
absorption and refraction (|Draind Il998t ) . Addition- 
ally, H I observations would reveal any associated neu- 
tral structures, Faraday rotation measurements may 
constrain the magnetic field strength within the struc- 
tures, and high resolution observations could confirm 
predicted changes in the apparent structure of the source. 
While there have been previous efforts in conducting 
such o bservations (|Clegg et al.lll996l : lLazio et al.ll2000l , 
bOOlal ). only one source has been observed in such a 
manner (B1741— 038), and, for some of the measure- 
ments, only a single observation during the actual ESE 
could be obtained. An all-sky survey could detect mul- 
tiple ESEs during the course of a year, allowing am- 
ple opportunity to trigger follow-up observations with 
other telescopes and at other wavelengths. 

2.2.3 Interstellar Scintillation of AGN 

Observations of decimetre v ariability in extragalactic 
sources (e.g.. lHunsteadiri972l ^ led quickly to the recog- 
nition that the sources needed to have extremely com- 
pact compone nts, be affected by interstellar scintilla- 
tion , or both l|Armstrong et aP Il977l : ISpangler et al.l 
11989 ). The need for extremely compact components, 
apparently violating the so-called inverse Compton catas- 
trophe, beca me even more a cute with the discovery of 
"flickering" (|Heeschenlll984l ). and later intra-day vari- 
ability. 

The MASIV 5GHz VLA S urvev (iLovell et al.ll2003l. 

2008) and other observations ziora-Chudcze r et al.l 
1997; Dcimctt-Th orpe fc de B ruvn 2002; Bigna ll et all 
20031 . 12OO6I ) have established that the fast — intra-day 
and inter-day — flux density variations at centimetre 
wavelengths exhibited by a subset of compact quasars 
are predominantly caused by ISS rather than intrin- 
sic variability. Intra-day variability refers to sources 
whose characteristic timescale is less than 24 hours, 
while inter-day variability refers to sources whose timescale 
exceeds this. 

Scintillating AGN are of astrophysical interest be- 
cause the small angular sizes that they must possess 
in order to exhibit rapid ISS require brightness tem- 
peratures near, or possibly in excess of, the lO^'^ K in- 
verse Compton limit for incoherent synchrotron radi- 
ation. In several cases over the past decade, it ap- 
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peared that the Doppler boosting factors required in 
order to reconcile the apparent brightness tempera- 
ture with the inverse Compton limit were well in ex- 
cess of the typical values determined from Very Long 
Baseline interferometric observations of quasars. Es- 
timates of the Doppler boosting factors are subject 
to considerable uncertainty related to the distance of 
the scattering material from Earth, but many of the 
most rapid IDVs appear to require Doppler factors of 
a few tens (iMacguart et al.ll2000l: iRickett et al ll2002l : 
iBignall et al.ll2006l : iMacauart fc de Bruvnll20071 1. 

More recently, it appears that the fastest mani- 
festations of IDV are associated with sources whose 
lines of sight intersect more local, albeit much more 
inhomogeneous, patches of turbulence than previously 
supposed. Thus there has been a transfer of our ig- 
norance away from the physics of the AGN and onto 
the physics of interstellar turbulence. However, this 
issue is not yet conclusively resolved, as little is known 
about the origin (or physical plausibility) of the local, 
highly inhomogeneous scattering screens that need to 
be invoked to explain the observed scintillation prop- 
erties of these sources. 

The four-epoch MASIV survey (|Lovell et al.ll2003l . 

I2OO8I ') detected variability in 56% of the 482 sources ob- 
served but the characteristics of the variability changed 
from epoch to epoch. While able to illustrate this 
ISS intermittency, the MASIV survey itself, because 
of its limited duration and cadence, could not address 
many of the resulting questions, notably the extent to 
which the intermittency results from changes in the 
ISM properties or changes in the sources. In order 
to address the cause of the intermittency, a survey 
must have the following characteristics. It must ob- 
serve a large number of sources widely distributed on 
the sky so as to have robust statistics and to be able to 
search for correlations expected for ISS (e.g., a Galac- 
tic latitude dependence or with nearby structures in 
the ISM). The survey must measure both the total 
and linear polarization intensity of the sources — ISM- 
driven intermittency should result in the linear polar- 
ization scintillation tracking that of the total inten- 
sity, while source-driven intermittency should result in 
the linear polarization and total intensity scintillating 
nearly independently of each other. The survey must 
be conducted with close to daily cadence to capture 
the scintillation variations, and over a multi-year time 
scale in order to separate changes in source structure 
or scattering medium from changes associated with the 
annual cycle resulting from the Earth's orbital velocity. 

In addition to determining the cause of ISS inter- 
mittency and mapping the scattering properties of the 
local ISM, such a survey would address the existence of 
AGN emission above the inverse Compton brightness 
limit and potentially constrain the baryonic content of 
the IGM as a function of redshift. 

The cause of ISS intermittency: By continuously 
observing a large sample of sources over a long pe- 
riod, VAST will comprehensively determine the inter- 
mittency properties of scintillating AGN. It will dis- 
criminate between ISM-related and source-driven ef- 
fects via polarization variations, as discussed above. A 



complication arises because annual cycles in the vari- 
ability timescale (caused by the Earth's velocity vec- 
tor changing with respect to the scattering screens) 
need to be distinguished from true IDV intermittency. 
However, these two phenomena can be separated by 
making regular (daily) observations over the course of 
at least two years. The variations should be dissimilar 
between epochs separated by a year if the scintillation 
is intermittent. Intermittency, if caused by the ISM, 
points to the physical nature of the turbulent patches 
responsible for IDV. 

AGN emission above the inverse Compton bright- 
ness limit: In some cases fast scintillation may be at- 
tributable to the appearance of micro- arcsecond struc- 

ture within a source, associated with an outburst (|Macauart fc de'fi' 

120071 '). The rise and fade times of the manifestation 

of fast variations allow us to determine the longevity 

of micro- arcsecond components associated with IDV: 

The rate at which the AGN supplies energy to power 

the bright emission that causes the variability, and the 

rate at which energy losses cause its eventual decay. 

These rates can then be compared with specific AGN 

processes to identify the mechanism associated with 

the super-Compton radio emission. 

Baryonic content of the intergalactic medium: 

The MASIV surve y showed that few sources above 
z ~ 2 are variable (|Lovell et al.ll2008l ). If this appar- 
ent suppression of ISS as a function of source redshift 
is confirmed by VAST, there are a number of possible 
explanations: (1) intrinsic source evolution; (2) gravi- 
tational microlensing; or (3) scattering in the turbulent 
ionised IGM. If the last is at least partially responsible, 
then the angular broadening is caused by the cumula- 
tive effect of all the baryons in the ionised IGM. The 
sheer number of sources VAST can monitor will enable 
us to probe in detail the evolution of structure in the 
ionised IGM as a function of redshift. 

One potential complication with scintillation stud- 
ies is that the largest amplitude IDV is typically ob- 
served at frequencies near 5 GHz. At ASKAP frequen- 
cies, it is possible that the strength of scattering, even 
toward high Galactic latitudes, and intrinsic source 
diameters will r esult in large-amp litude scintillations 
being quenched (jDennison fc Condon, 1981, '). although 
lower amplitud e refractive scintillat ions may still be 
exhibited (e.g.. iDennison et al.lll987l '). 

2.3 Accretion and Magnetism 

Along with explosions and propagation effects, accre- 
tion and magnetic fields are frequently implicated in 
transient events. Accretion onto a compact object is a 
powerful energy source that drives both local sources 
such as X-ray binaries (XRBs) and microquasars as 
well as extragalactic blazars and AGN at cosmologi- 
cal distances. The accretion flow is often connected to 
an outflow in the shape of a jet, which is the source 
of the radio emission. The intimate link between the 
accretion disc and the jet is demonstrated by the fun- 
damental plane of black hole activity that describes a 
correlation between the mass of the black hole, the jet 
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power, and the ac cretion power (jMerloni et al ] |2003l : 
iFalcke et al.ll20d^ in active galactic nuclei. This re- 
lation can be extended to lighter black holes where 
it matches the correlation between the X-ray and ra- 
dio luminosities in XRB black hole ca ndidates (BHCs) 
that was found bv lGallo et aD (|2003l l. Radio emission 
is also observed from XRBs containing weakly mag- 
netised neutron stars, such as Z and atoll sources, al- 
though they are at least one orde r of magnitud e fainter 
in radio than the BHCs (e.g. lFendejT2006l '). and a 
weak radio jet has also been observed from the cata- 
clysmic variable SS C yg during a dwarf nova outburst 
(|Kording et al.|[200^ . There is still not a good theo- 
retical understanding of how accretion discs produce 
outflows, but most models assum e that the jet is ac- 
celer ated by a magnetic field (e.g. iBlandford fc Pavnd 
119821 '). VAST will aUow us to characterise the full 
range of transient jet phenomena across the fundamen- 
tal plane. 

Magnetic fields can also act in other ways. They 
can convert the rotational energy of their source into 
electromagnetic energy and the mechanical energy of 
an outflow, as happens for instance in a radio pulsar, or 
they can release their own energy in the form of flares 
on a star or the radio emission from Jupiter and (pos- 
sibly) Jupiter-like extrasolar planets. With VAST we 
will be able to reveal the unifying physical principles 
that underlie such diverse behaviour. 



2.3.1 Intrinsic Variability in Active Galac- 
tic Nuclei 



As discussed in Section lTT^lFiedler et al.l (llQSTaT l ob- 
served a sample of 33 AGN, with flux densities of a 
few Janskys, every day for more than 6 years. All 
sources were variable over this time span with typi- 
cally smooth variations in flux with quasi-periods of 
months to years. VAST will produce light curves and 
detect this sort of variability for more than 10* sources 
over the whole sky. The Fermi gamma-ray satellite ob- 
serves the whole sky every 3 hours and has detected 
blazar variability on timescales as short as 6 hours 
(jAbdo et al.ll2010l ). Simultaneous measurements of the 
gamma-ray and radio variability in blazars will provide 
us with a powerful tool for understanding the physics 
of the central engines in AGN. 

Beyond single black holes in AGN, a number of bi- 
nary supermass ive bla ck holes (SMBHs) arc know n or 
suspected (e.g., iKomossa 2006; Rodriguez ct al. 20061 : 
ISmith et al1l2010l : lBurke-Spolaoi)r2010l '). The evidence 
for these varies, but it is clear that with the mergers of 
galaxies the central SMBHs could also merge. The act 
of inspiralling may lead to radio transients accessible 
to ASKAP. This is because the SMBHs will likely move 
through a static magnetic field provided by an accre- 
tion disk, and this should drive electr o magnetic jets 
(iPalenzuela et al.|[2O10l : iLvutikovl [2OI1I : iMoesta et aD 
I2OI2I ). The detailed emission processes are unclear, 
but it is plausible that there could be a radio signature 
accessible t o ASKAP before (lO'Shaughnessv et al.ll201ll ') 
or during (|Kapl an et al. I I2OIII ) the merger. Numer- 
ous other ele ctromagnetic mech anisms have also been 
proposed; see lSchnittmanI (120111 ) for a detailed discus- 



sion. Combining the radio observations with observa- 
tions at other wavelengths or even with gravitational 
waves would significant ly improve our un derstanding 
of galaxy mergers (e.g.. ISesana et al]|201ll l. 

2.3.2 Tidal Disruption Events 

Tidal disruption events (TDEs) occur when a star wan- 
ders too close to a black hole at the centre of a galaxy 
or star cluster. The black hole tidal forces tear the 
star apart, but the debris continue on ballistic orbits 
around the star and are accreted as they approach pe- 
riastron on the following orbit. The resulting sudden 
increase in the black hole accretion rate may be di- 
rectly detected at optical and X-ray frequencies. After 
the in itial flare up the luminosity decays as t~^^^ (jReesI 
Il990l l. Several candidates of such TDEs have been de- 
tected serendip itously through all - sky X-ray and ultra- 
viole t surveys (i Grupc ct al. 199^; Komossa & Grcine: 
19991; IKomossa fc B ade 1999; Esquci ct al..,2007,.,20o' ^ 



Gczari c t al.ll2006l. 12008. 200 9l;lLin et al.l201]l ; lGezari et al.l 
20121 : iCenko et al.ll2012ai : iSaxton et al.lf2012r i . 



Crucially, black hole accretion is often accompa- 
nied by relativistic jets, with a mechanical power often 
comparable to the radiative luminosity. The jet can be 
a source of beamed radio emission and hard X-rays in 
a way similar to a blazar. To a hard X-ray telescope 
such as the Burst Alert Telescope (BAT) on the Swift 
it will initially resemble a GRB, but it will be distin- 
guished by its much longer duration. Swift J1644+57 
is the first observed example of such a TDE candi- 
date, at redshift z — 0.35. It was discovered as a lu- 
minous X-ray transient and was observed across the 



^^^^^^Q^^ — ^ ^^^^^^^^^ 

ysis and modeling of the event ( 


Bloom et all 1201 ll: 


Levan et al.l 12011': 'Zaudoror et al.l 


201 ll: iBerger et all 


2OI2I: IWiersema ct al. 2012,). The mass of the cen- 



tral black hole responsible for this event was derived 
to be < lO'^M© and the ejecta Lorentz factor > 10. 
More recently another TDE candidate has also been 
detected by Swift in hard X-rays, Swift J2058.4+0516 
(|Cenko et al.ll2012bl ). 

The main difference between AGN and TDE jets is 
that the former propagate to Mpc distances, and cre- 
ate characteristic cocoon/lobe structures where their 
kinetic power is dissipated. By contrast, TDE jets are 
impulsive ejections that produce a forward sh ock /reverse 
shock structure as they interact with the ISM l|Sari fc PiranI 
1 19951 '). As there is no pre-existing cavity around the 
previously quiescent black hole, TDE jets have to prop- 
agate through a much denser environment, and de- 
celerate to sub-relativistic speed at distances < 1 pc. 
Most of their kinetic power is expected to be dis sipated 
in the reverse shock (|Giannios fc Metzgeill201ll ). Typ- 
ical TDE radio flares from the disruption of a solar- 
mass star by a 10^ solar mass black hole are expected 
to peak in brightness after ~ 1 year, reaching a flux 
density of ~ 2 mjy at 1.4 GHz, for a distance of 1 
Gpc corresponding to 2 = 0.2 (not including the ef- 
fect of relativistic boosting for jets oriented along our 
line of sight). At peak brightness, the expected ra- 
dio spectrum is inverted (self-absorbed) with F,^ ~ 
2 GHz, and F^ ~ above that 
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(|Giannios fc MetzgedbOHl ). 

Radio afterglows of TDEs would be a powerful 
tool to understand jet physics and source geometry, 
through measurement of the flux evolution and po- 
larization. The expected degree of polarization for 
these events is low if they are collimated, similar to 
GRB afterglows. Recent infra-red and radio polariza- 
tion measurements have been found to be consistent 
with theoretical ex pectations (|Wiersema et al. 1 120121 : 
iMetzger et"al]|2012r i. 

Radio detections of TDE jets hold the key to other 
unso lved astrophysical problems (e.g. iKrolik fc PiranI 
I2OI2I I. It is well known that not all accreting black 
holes have jets, but it is not clear why. Jets are formed 
only at certain ranges of accretion rate, but other pa- 
rameters such as the black hole spin may also be im- 
portant. Finding what fraction and what type of X-ray 
detected TDEs have radio jets will help address this 
question, and will shed light on the evolution of nuclear 
black holes in galaxies. Furthermore, the timescale for 
the formation of a radio-bright TDE jet will help us 
understand how quickly a collimated magnetic field is 
generated inside the transient accretion disk. 

The predicted brightness of TDEs at gigahertz fre- 
quencies suggests that VAST should be able to detect 
about a dozen events per year. A systematic survey 
of the radio sky would be able to constrain the rate of 
TDEs and reduce the coUimation uncertainty signifi- 
cantly. In particular, the discovery of an entirely new 
population of orphan TDEs would allow us to place 
solid constraints on the occurrence rate of this phe- 
nomenon and therefore test the theoretic al predictio n 
that we should see ~10~'* yr~^ per galaxy (|ReeJl990l '). 

2.3.3 Accreting Neutron Stars, Black Holes 
and Microquasars 

Exploring the connection between accretion disks and 
jets in X-ray binaries (XR Bs) requires c orrelated ra- 
dio and X-ray observations (|Fendei1l200^ ') . The VAST 
survey will reveal and subsequently measure the radio 
emission from a large number of XRBs and allow us to 
search for new radio outbursts. This, coupled with the 
recent availability of large field of view instruments like 
Swift-'BAT and Fermi-GBM, and with the early results 
(from 2013) of eROSITA's all sky survey, will provide 
a large sample of XRBs with active jets, and will al- 
low us to do statistical studies of their time variability 
properties and X-ray/radio correlation. 

VAST will yield information on the population of 
XRBs showing quasi-steady state radio jets, which are 
mildly relativistic, and are associated with non-thermal, 
radiatively-inefiicient accretion states ([Fender et al.|[20M : 
iMcClintock fc RemillardI 120061 ). For these persistent 
jets, fundamental properties such baryonic content, en- 
ergy budget, emission mechanism, and the role played 
by the s pin of the compact object are not well con- 
strain ed (|Markoff et al.ll2005l : lNaravanll99^ : lYuan et al] 
l2005l l. and regular radio monitoring will allow us to 
investigate them. In addition, we will determine the 
duty cycle of the radio-loud phases, and test whether 
radio state transitions (changes in the kinetic power) 
are always associated with X-ray transitions (changes 



in the radiative power). We will also investigate what 
physical mechanism causes major radio flares and the 
ejection of highly relativistic, optically-thin bullets, of- 
ten associated to the state transition from the hard/non- 
thermal to soft /thermal state. 

In addition to stellar mass and supermassive black 
holes, it is predicted that VAST will detect jets from 
intermediate mass black holes (IMBHs) with masses 
between ~100 ~ 100,000 Mq. IMBHs may provide 
a pathwa y for the formation o f supermassive black 
holes (e.g. lEbisuzaki et al.ll200ll ) and also have impor- 
tant connotations for other areas of astrop hysics (e.g. 
iFornasa fc Bertone|[200i : IWang et al.ll2010l ). but until 
recently their existence was highly disputed. The dis- 
covery of the object known as HLX-1 located in the 
galaxy ESO 243-49, with an X-ray luminosity ~1000 
times greater than the Eddington limit for a 10 M0 
black hole, currently pro vides very stro ng evidence for 
the existence of IMBHs (|Farrell et al.|[2009. ). with in- 
dependent studies constraining the mass to between 
9,000 M(7) < Mbh < 90,000 Mr (iDavis et al]|201ll: 
IServillat et al.l201ll : iGodet et al.ll2012l : IWebb et al.l2012l ) . 
Recently, transient radio emission was detected from 
HLX-1 following a transition from the low/hard to 
high/soft spectral states, consiste nt with the ejectio n 
of relativistic jets from an IMBH (|Webb et al.ll2012r ). 

Objects such as as HLXl may be very rare (only 
one known within 100 Mpc), but there a re hundreds of 
ultraluminous X-ray sources (ULXs; see lFeng fc Sorial 
I2OIII for a recent review) in the local Universe with 
luminosities exceeding the classical Eddington limit 
of stellar-mass black holes. Most of them are likely 
to be explained through other mechanisms such as 
anisotropic emission and/or super-Eddington accretion, 
but we cannot rule out that a few of them may be 
IMBHs. In addition, there may be hundreds of low- 
state or quiescent IMBHs lurking in the core of glob- 
ular clusters, in the nuclei of dwarf galaxies, or inside 
recently merged galaxies that may occasionally become 
active. The detection of jet emission from these objects 
would allow us to infer the black hole mass, as the ra- 
dio and X-r ay luminositie s are correlated as a function 
of mass (e.g. lFendej|2006l ) . Steady jets would be unde- 
tectable for all but the closest of these objects. How- 
ever, transient radio flares have been observed from 
XRBs associated with jet ejection events with radio lu- 
minosities a factor of 10 - 100 ti mes brighter than th e 
non-fiaring radio emission (e.g. iKording et alll2006l ). 
Such fiares should be detectable by VAST out to dis- 
tances of -^5 - 15 Mpc from 1,000 M0 black holes, and 
out to ~23 - 72 Mpc from 10,000 M0 black holes (Far- 
rell, in preparation). Thus we will be able to survey the 
local Universe for flaring IMBHs, placing constraints 
on the number density and therefore testing their va- 
lidity as a pathway for supermassive black hole forma- 
tion. Conversely, if ULXs are super-Eddington accre- 
tors, detection of radio emission will prove that even 
at those extreme accretion rates, a substantial fraction 
of accretion power can come out as kinetic power, and 
will improve our understanding of black hole feedback 
at high accretion rates. 
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2.3.4 Neutron Star Mergers 

The merger of two compact objects (neutron stars) is 
expected to generate strong gravitational waves (' i)2.1.3|) 
but they wiU also reveal themselves through electro- 
magnetic radiation. The relativistic, beamed outflows 
detected from short-hard G RBs are thought to origi- 
nate from such mergers fe.g lFong et alll2012l ). VAST 
is not expected to fortuitously observe such beamed 
events. However, merger simulations also show that 
substantial mass can be ejected from these systems 
quasi-isotropically and at sub- or trans-rel ativistic ve- 
lociti es, and with energies of 10^^-10^° erg (|Piran et al.l 
120121 '). 

Consider a mildly relativistic blast wave (/3o ~ 0.8) 
moving out with an energy ~ 10*^ erg into the sur- 
rounding medium. The shock- wave sweeps the sur- 
rounding medium and heats it to relativistic tempera- 
tures by converting bulk kinetic energy of the incom- 
ing material into thermal energy of the shocked ma- 
terial. The electrons in the shock-heated plasma are 
accelerated in the post-shock magnetic field to radi- 
ate synchrotron radiation. This would eventually be 
optically thin in the gigahertz regime, peaking on a 
timescale of about a couple of months when the blast 
wave starts to decelerate. The optimum waveband for 
detecting these transients, which are also supposed to 
be prime sources of gravitational wa,ves, i s around gi- 
gahertz frequencies (|Nakar fc Piranlboilh . A merger 
within 300 Mpc would shine as a radio transient of 
about 1 mjy brightness at ~1 GHz (Kamble & Ka- 
plan in preparation) which would be detectable to the 
VAST and, potentially, also to Advanced-LIGO. Both 
the peak fiux density and the time to peak are un- 
certain since we lack accurate estimates of the ejecta 
velocity and energy, and we expect the circumburst 
density to vary from 1 cm"'' to 10~* cm~^. 

Using the 'best-bet' neutron star merger detec- 
tion r ates predicted for Advanced LIGO (|Abadie et al.l 
120101 ) ■ we estimate that up to about half a dozen of 
these transients should be detectable by ASKAP at 
any given time, although there are significant uncer- 
tainties involved. iPiran et al.l (|2012l ) independently es- 
timate ASKAP's sensitivity to these events and obtain 
similar results. These are comparable to the rate es- 
timated from observations of short GRBs which have 



been pr oposed to originate i n sim ilar mergers (jCoward et 
I2OI2I '). iMetzger fc Bergeil (|2012l ') emphasise the chal- 
lenges related to uncertain GW event localisations, 
which are significant. As discussed in i)2.1.3l the large 
instantaneous sky coverage of ASKAP may be a sig- 
nificant advantage in this situation. 



2.3.5 Flaring and Pulsing Neutron Stars 

Of the known neutron stars, the vast majority are pul- 
sars that have been detected through their pulsed radio 
emission. VAST will not be sensitive to radio pulsa- 
tions on timescales of < 1 second, but will instead 
probe neutron star populations that are not selected 
in pulsar surveys. We will be sensitive to intermittent 
or deeply nulling radio pulsars as well as radio pulsars 
with large variations in their average pulsed flux. 



Magnetars will be detectable through their rare 
flaring events such as the giant flare an d expanding ra- 
dio nebula produced by SG R 1806—20 (jGaensler et al.l 
I2OO5I : lUameron et al.ll2005h . In addition, some magne- 
tars are known to turn on as radio sources. For exam- 
ple, XTE J1810— 1 97 was initially dete cted as a tran- 
sient radio source (|Halpern et al.l 120051) before i ts ra- 
dio pulsations were identifled (|Camilo et ai]|2006l ). We 
may also be able to detect bright intermittent pulses 
from the enigmatic rota ting radio transients (RRATs; 
iMcLaughlin et al.l2006l ). VAST will thus reveal a broader 
census of the Galactic neutron star population than 
a standard pulsation search, with subsequent impli- 
cations for the overall supernova event rate in our 
Galaxy. 

The radio source FIRST J1023-I-0038 presents an 
interesting case study. Initially it was classified as a 
Galactic cataclysmic variable (CV) ([Bond et al.ll2003 ) 
because of its optical spectrum and variability, but its 
radio flaring behavior appeared to be anomalous. The 

classiflcation was called into question because of the 

chan ges in the spectrum over time ('e.g.. lThorstensen fc Armstrong 
I2OO5I ) and it was eventually re-identified as a transi- 
tional object where a recycled radio pu lsar had turned 
on af t er a lo w mass X-ray binary phase ([Archibald et al.l 
l2009l . [2010l ). The actual transition event had an ex- 
tremely short duration, with the accretion disk ap- 
parently dissipating between 2001 December and 2002 
May. The radio pulsar was probably absent in 1998 
(upper limits of 1.8 and 3.4 mJy at 1.4 GHz) but de- 
tected in a survey in mid-2007 with a flux density ~14 
mJy at 1.6 GHz. 

2.3.6 Classical novae 

Novae occur in binaries (CVs) that consist of a white 
dwarf that is accreting matter from its companion star. 
A nova is a thermonuclear explosion in the layer of ac- 
creted matter on the surface of a white dwarf. A shell 
of matter is ejected during the explosion and this shell 
might produce radio emission. The radio emission can 
persist for a few years after the explosion, and in the 
case of FH Ser it reached ~10 mJy at 2.7 GHz about 
one year after the outburst. In the late stages when 
the ejecta is optically thin, radio emission is generated 

1 by the entire ejecta and can therefor e be used to es- 

■^timate the total ejecta mass (see e.g., lHiellming|[l990l 
for a review). 

2.3.7 Flare stars, pre-main sequence stars, 
and rotation-driven stellar activity 

There are several types of stars that display large flares 
that can be observed in the radio band. The most com- 
mon are cool dwarf stars of spectral type M or later 
including brown dwarfs. The flares have typical flux 
densities above ~ 1 mJy, but there is also quiescent 
emission at the sub-m jy level between the flares (e.g. 
iPestalozzi et al.ll2000l ). For these types of systems it 
has been argued that radio observations are the best 
way to ob tain informat ion on the stellar magnetic field 
strength (|Bergei]|2006l ). The number of fiaring events 
that are potentially detectable is unclear, but iBerged 
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(|2006l l obtained a detection rate of ~ 10% for objects 
of later spectral class than M7. They also find that the 
radio to X-ray luminosity ratio increases over at least 
three orders of magnitude for these objects compared 
to earlier spectral classes. VAST will undertake a cen- 
sus of such stars in the local (10 pc) neighbourhood. 

In AM Her-like systems (cataclysmic variables in 
which the white dwarf has a surface magnetic field of 
~ 10^ - 10" T) the M dwarfs are forced to co-rotate 
with the binary, which has an orbital period of a few 
hours. The dwarf star then becomes more active; flares 
in AM Her itself h ave been observed to be as bright as 
10 mJy at 5 GHz (|Dulk et al.lll983l l. 

Another group of stars t hat can produce ra dio fiares 
are the RS CVn-systems (|Gunn et all llQQ^ L These 
are close detached binaries consisting of two main se- 
quence or evolved stars of spectral type F or later, that 
are tidally locked. The rapid spin of the stars leads to 
strong magnetic activity, which is manifested as large 
spots on the stars and also as fiaring activity that can 
be observed both in X-rays and at radio wavelengths. 

Flare-like outbursts are also observed in pre-main 
sequence stars, but it is not clear whether the physical 
mechanism is really that of a stellar flare or related 
instead to the accretion flow onto the star. Such out- 
bursts have been observed in both naked T Tauri stars, 
which do not show signs of extended accretion discs, 
and in systems with an infrared excess, whic h indicates 
the presence of a disc (|Osten fc Wolkll2009l ). 

A magnetic field can convert the rotational energy 
of a star into an outflow and electromagnetic radiation. 
The most well-known examples of this are the classical 
radio pulsars, but a similar phenomenon can arise in 
less extreme stars. Radio observations have detected 
two fully circularly polarised radio pulses per spin pe- 
riod from the Ap star CU Vir, which is also spinning 
down in a way similar to a pulsar (e.g . iTriglio et al.l 
120081 : iTrigilio et ai]|201ll : iLo et all [2013 ). 

Another example of activity driven by rotation is 
the white dwarf in the cataclysmic variable AE Aqr. 
The spin period of the white dwarf is 33 s, so rapid 
that the white dwarf does not accrete the matter that 
streams over from the companion K dwarf. Instead, 
the matter is ejected from the system. AE Aqr can be 
observed over a wide range of frequencies from X-rays 
to radio, and pulsed emission is observed at optical and 
X-rays, but it has not yet been detected at radio wave- 
lengths (iDe Jaeer et al]| 1994 [Bookbinder fc LambI 19871 : 
iBastian et al.ll 19881 ). 



3 Exploration of the Unknown 

The most interesting sources detected by VAST may 
be those we currently know nothing about. A sensi- 
tive wide-area blind survey is ideal for detecting such 
unknown source classes, but one of the challenges in 
designing the VAST survey is ensuring that the sur- 
vey parameters do not unduly bias us against unknown 
source classes that we would otherwise be sensitive to. 
For example, daily monitoring of a field for 5 minutes 
per day makes it less likely that we will detect an ob- 
ject that brightens on hour-long timescales and fades 




ft are 



flare 



-8-6-4-2 Q 2 4 6 8 10 
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Figure 2: The param eter space fo r radio tran- 
sients, adapted from [Cordes et al.[ ( 2004) . A 
quantity equivalent to absolute luminosity (ob- 
served flux density S multiplied by the square 
of the distance D^) is plotted against the di- 
mensionless product of the emission frequency v 
and the transient duration or pulse width W. In 
the Rayleigh- Jeans approximation, these quanti- 
ties are directly proportional and related to the 
brightness temperature, as indicated by the diago- 
nal lines, with T ~ 10^^ K marking the maximum 
brightness temperature of incoherent processes. 
Examples of known radio transient sources are in- 
dicated, including giant pulses an d "nanogiant" 
pulse s from the Crab pulsar (iHankins et al.[ 
l2003fl radio puls es from XTE J1810-197 
([Camilo et al.l [20061 ) and other pulsars from the 
Austral ia Telescope National F acility pulsar cat- 
alogue (jManchester et al.][2005[): the niicroq uasar 
GRS 1915-H05 ([Mirabel fc R odn'gued [19941): ra - 



dio flares from V4641 Sgr ( Hiellming et al 
the brown dwarf LP 944-20 (Bereeretal 



2000 ). 



2001), 



rge 

and the m agnetar SGR 1806—20 (* Gaensler et al.[ 
[2005[ : [Cameron et al.. ^005) ; the Galactic centre 



radio transient J1745— 3009 (jHvman et al.l l 
puls es from the ultrac ool dwarf TVLM 513—46546 
( Hal linan et al.[ [2003) ; as well as radio emission 
from the Sun and Jupiter. Red lines indicate the 
expected sensitivity of ASKAP to sources at dis- 
tances of 10 pc, 1 kpc and 1 Mpc. 



within several hours. 

Large-scale blind surveys for radio tran sients have 
been performed by a number of groups (e.g. .lAmv et al.l 
[19891 : iKatz et al.[ [20031 : [Matsumura et alj [2007[ ) over 
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the last few decades. These surveys have been severely 
limited in their sensitivity, sky coverage, and cadence. 
In spite of this, they have revealed that the radio sky 
contains many transient objects, the identities of some 
of which remain mysterious. 

For example. I Hvman et al.l (|2005l ) discovered a burst- 
ing transient towards the Galactic Centre which lasted 
for only a few minutes with a flux density in excess of 
2 Jy at 330 MHz. Subsequent observations showed 
that the bursts repeated multiple times with a period 
of ~ 77 minutes and a burst length of ~ 10 minutes. 
The identification of this source re mains unclear (e.g. 
iKaplan et alll2008l : IRov et al.ll2010l '). al though similar- 
ities with bursts from ultracool dwarfs (|Hallinan et al.l 
120071 ') are intruiging. Another example is the well- 
known Lorimer burst, discovered in a n arc hival search 
of Parkes pulsar data. iLorimer et al.l (|2007f ) discovered 
a single millisecond pulse of extragalactic origin with 
the astonishing peak flux density of 30 Jy at 1.4 GHz 
(also sec Kcanc ct al. 201ll for a second example). Al- 
though the nature of this source has been questioned, 
and recent work suggests t he burst may be terrestrial 
(|Burke-Spolaor et al.|[201ll ). these new discoveries il- 
lustrate that the parameter space for radio variability 
has not been fully explored yet (Figure [21). 

In this section we summarise existing blind and 
archival surveys for radio transients, focusing on 'slow' 
transients with variability on timescales of seconds or 
greater. 

3.1 Transient Rates from Blind Sur- 
veys 



view Q. We can rearrange Equation [2] in terms of p 
for the purpose of placing theoretical constraints on 
planned surveys: 



iBower et al.l (|2007l ) presented the results from an archival 
survey of VLA data. Since the publication of this work 
there have been a number of blind radio transient sur- 
veys conducted on both archival data, and from the 
flrst datasets produced by next generation radio tele- 
scopes. Most of these have had the dual aim of (a) 
searching for individual objects of interest and trying 
to identify them; and (b) characterising the population 
statistics of radio transients in order to prepare for fu- 
ture surveys. Although the reported transient rates 
from iBower et al.l (|2007l ) have recently been rev ised, 
after a reanalysis of the data bv lFrail et al.l (|2012l '). the 
work established a framework for comparison of radio 
transient surveys. 

I Bower et al.l (|2007l ) proposed a metric for compar- 
ing the results from blind surveys for radio transients: 
the two-epoch equivalent snapshot rate, which is equiv- 
alent to the transient surface density. Figure [3] shows 
an updated version of their plot, incorporating a sub- 
set of the surveys listed in Table O and discussed in 
the rest of this section. For surveys that made no de- 
tection, an upper limit is calculated using a Poisson 
distribution: 

P(n) = e-""^ (2) 

where P(n) is the confidence level (i.e. P{n) = 0.05 for 
95% confidence), p is the areal density of transients in a 
two-epoch survey, equivalent to the snapshot transient 
rate, and A is the equivalent solid angle, calculated by 
multiplying the number of images Nt by the field of 



InP(n) 
nx Nt' 



(3) 



We can improve Equation[3]by taking into account 
not only the total number of observations in a given 
survey Nt, but the intra-observation cadence r, out 
of total observing time per observation T. Doing this 
results in: 

p = 'iim . (4) 

For example, on the smallest timescale, ASKAP im- 
ages will be produced every r = 5 seconds; the total 
observing time per observation (in one scenario) could 
be T = 1 hour and the field could be observed daily for 
two years {Nt = 730). A five second ASKAP image 
will have an RMS of 1.4 mjy beam~^. If subsequent 
data or images are stacked together (either via the uv 
or image plane) the RMS will improve to: 



1.4 mJy beam ' 

^Nt X (T/rf 



(5) 



Equations [4] and [5] have been used to generate sur- 
vey predictions from the VAST survey parameters de- 
scribed in Section|4l and these predictions are shown in 
Figure [3] We also include the pre dicted rates of a s ub- 
set of known radio transients (see iFrail et aL l l2012l . for 
futher discussion). One limitation of the plot in Fig- 
ure[3]is that it collapses a number of dimensions into a 
single two-dimensional plot. For example, it does not 
represent the cadence of the surveys, or any trends in 
transient detection rate with observing frequency. 
Surveys typically fall into one of two cat egories: 



many repeated observati ons of a single field fe.g. lCarilli et al] 
12003 : [B ower et all 120071'), or a large area, few-epoch 
surveys Ce.g. lLevinson et al.l2002l : iGal-Yam et~al]|2006l ') . 

The two-epoch equivalent snapshot rate is a valid met- 
ric for quantising the effectiveness of a survey, pro- 
vided the the timescale of transient behaviour is less 
than the cadence of the observations but longer than 
a single integration. In this case, a survey with many 
epochs is equivalent to a larger area survey with two 
epochs. Further survey metrics have been developed to 
account for observation cadence and source evolution 
timescales as well (e.g.. ICordesll200^ ). 

3.2 Archival Surveys 

The first large radio archival survey of variability and 
transient behaviour was conducted bv lLevinson et al\ 
(|2002l l who covered a large fraction of the sky with 
two epochs by co mparing the NVSS (|Condon et al.l 
1 19981 ') and FIRST (| White et al.lll997l ) surveys, both 
at 1.4 GHz. Various corrections were made to ac- 
count for the different footprint of each survey and 
their differing resolutions (45" for NVSS compared to 
5" for FIRS T). Nine possi b le rad io transients were 
reported by iLevinson et al.l (|2002|). but follow-up o f 
these sources with the VLA bv lGal Yam et all (|2006l ) 
ruled out five of these as false candidates and two as 



Table 3: Summary of snapshot rates for transient and variables radio sources reported in the literature. The results are organised according to upper 
limits based on non-detections (top section); transient detections (middle section); and detections of highly variable radio sources (bottom section). 



Survey/Reference 




^Max 


Rate 


Timcscale 


Frequency 


Epochs 




(mJy) 


(mJy) 


(dcg-2) 




(GHz) 


(N) 


Bower et al. (2007) 


>0.09 


— 


<6 


1 year 


4.8 & 8.4 


17 


Bower et al. f2007) & Frail et al. f2012)t 


0.2 


— 


<3 


2 months 


4.8 & 8.4 


96 


Bower et al. f2007) & Frail et al. f2012)t 


>0.37 


— 


< 6 X 10-1 


20 mins - 7 days 


4.8 & 8.4 


944 


PiGSS-I /Bower et al. f2010)('A)* 


>1 


— 


<1 


1 month 


3.1 


75 


PiGSS-II /Bower et al. ('2011VB)* 


> 5 


— 


<0.18 


1 month 


3.1 


78 


FIRST-NVSS /Gal- Yam ct al. (2006) 


>6 


— 


<1.5xl0"3 


days to months 


1.4 


2 


Beh et al. (2011) 


>8 (8cr) 


— 


<0.032 


4.3 - 45.3 days 


1.4, 4.8 & 8.4 


5037 


PiGSS-I/Bower et al. (2010) (A)* 


>10 


— 


<0.3 


1 month 


3.1 


75 


PiGSS-II /Bower et al. (2011)(B)* 


> 15 




<0.025 


1 day 


3.1 


78 


AlAlb - 1/Urott et al. (2010) 


>40 




<0.004 


81 days - 15 years 


1.4 


12 


r5ower iz Daui izUiillA) 


> ( U 




<oX iU 


1 day 


1 A 


iooz 


ATATS - II /Croft et al. (2011) 


>350 




<6xlO-4 


minutes to days 


1.4 


12 


Bower & Saul (2011)(B)* 


>3000 




>9xl0"'^ 


1 day 


1.4 


1852 


Lazio et al. (2010) 


>2.5xl06(5cr) 




<9.5xl0-** 


5 mins 


0.0738 


-1272 


JA12'^ Jaeeer et al. (2012) 


2.1 




<0.12 


1 day to 3 months 


0.325 


6 


Ofek'^/Ofek et al. (2011) 


2.4 




0.039 


1 day-2 years 


4.9 


16 


MOST"^ /Bannister et al. (2011a)* 


14(5cr) 


6500 


1.5 X 10-3 


days-years 


0.843 


3011 


M09'^/Matsumura et al. (2009) 


1000 


4300 


3x10-3 


- 1 day 


1.4 




C03^/CariUi et al. (2003) 


>0.1 




<18 


19 days and 17 months 


1.4 


5 


F03^ /Frail et al. (2003) 


>0.25 




5.8 


~ 1 day 


5 & 8.5 




Becker^ /Becker et al. (2010) 


>1 




1.6 


~ 15 years 


4.8 


3 


GT86^/Gretrorv & Tavlor (1986) 


4.6 




1x10-3 


1 day - 5 years 


4.9 


16 


MOST^ /Bannister et al. (2011a) 


>14 




0.268 


days-years 


0.843 


3011 



I^Note that althou gh iBower et al.l (l2007t ) report 10 detections, the reanalysis bv lFrail et al.l ( 2012| ) results in only 4 (uncertain) detections. The snapshot 

rate calculated by Frail et al.l (20121) ass umes no detections, and w e adopt this as the more conservative position for our analysis. 

iBower et al.l ( 20101 ) . Bower et al.l (|201l[ ) and lBower fc Saul ( 201lh each state two different rates depending on flux density and timescale; we quote 

these separatel y as (A) and (B). 

^See also iBannister et al.l ( 2011b[ ). 
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Figure 3: Log 2-epoch snapshot transient rate (deg~^) against log of the flux density (Jy) for surveys 
that report detections of transient and variables (thick lines) , we also include a selection of surveys that 
report upper limits (thin lines with arrows) — see Table [3] for survey acronyms. The surveys are coloured 
according to frequency: black < 1 GHz, blue = 1 —4.8 GHz and red > 4.8 GHz (see Table|3]for details). No 
corrections have been made for spectral index effects. Surveys labelled superscript 'V denote detections of 
highly variable radio sources; those with superscript 'T' denote detections of transient type sources. The 
VAST survey predictions arc indicated and organised by sub-survey. A description of each survey is given 
in Section m In each case the vertical segment denotes the RMS that can be achieved per observation. 
The horizontal segment indicates the upper limit that would be set if no transients or variables were 
detected in the entire survey. The joining line (between the horizontal and vertical) indicates the upper 
limits that could be placed by combining data from multiple five second snapshots up to the total available 
integration time per observation (see Table [21) ■ See Section [XT1 for further discussion. 



non-variable sources. The effective survey area was 
5990 deg^, hence with two detections the overall tran- 
sient snapshot rate was p = 1.5 X 10-3 deg-^ 

An other large area survey was conducted bv lBannister 
(|2011al 'l. using archival data fr om the Mo longlo Ob- 
servatory Synthesis Telescope (|Millsl IT"98lh which op- 
erates at 843 MHz. This dataset had the advantage 
of a common resolution and comparable sensitivity in 
each epoch, but like the NVSS-FIRST comparison had 
a ran ge of cadences b e tween days and years for each 
field. iBannister et all (|2011al lbh detected 2 transient 
sources from a survey with a two-epoch equivalent 



area of 2776 deg^, giving an overall snapshot rate of 
p = 1.5 X 10"'^ deg-^ for sources above 14 mjy. 

An alternative approach has been to use archival 

etcffllj brator fields which have many repeate d observa- 
tions but a relatively small field of view. iBell et al.l 
(20lT|) analysed 5037 observations of VLA calibrator 
fields at three frequencies (1.4, 4.8 and 8.4 GHz). The 
observations ran from 1984 to 2008, with typical gaps 
between observations of between 4.3 and 45.3 days. 
No radio transients were detected and they therefore 
placed an upper limit of p < 0.032 deg"^ for sources 
brighter than 8.0 mJy between 1.4 and 8.4 GHz. 
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Likewise iBower fc Saul l|201lh detected no tran- 
sients in their analysis of archival observations of 3C 286 
from the VLA. Their dataset spanned 23 years and in- 
cluded 1852 epochs at 1.4 GHz. They derive upper 
limits of p < 3 X 10""^ deg~^ at 70 mjy for timescales 
of ~ 1 day and p < 9 x 10~* deg~^ at 3 Jy for 
timescale s of ~ 1 m i nute. These results, combined 
with the ICroft et al.l (|2010t ) results discussed in Sec- 
tion 13.31 below, currently provide the best limits on 
the number of slow, bright radio transients. 

3.3 Gigahertz Surveys 

There have been a number of blind imaging surveys 
specifically designed to search for transient and highly 
variable sources. Until recently these have generally 
been small area surveys with many epochs, due to the 
limited field of view of the telescopes used. Early sur- 
veys were motivated by searching for radi o counter- 
arts to GRBs and supernovae. For example. iFrail et al.l 
19941 ) observed 15 epochs of a field centred on the lo- 
cation of GRB 940301 with the Dominion Radio Astro- 
physical Observatory Synthesis Telescope at 1.4 GHz. 
These observations were taken between 3 days and 99 
days after the initial event and made no detections, 
which resulted in an upper limit of 3.5 mJy on any 
ti me- variable radio sou rces . 

ICarilli et all (|2003l l observed the Lockman Hole re- 
gion to a sensitivity limit of 100/iJy at 1.4 GHz with 
five epochs over a period of 17 months. They found 
that 2% of the sources within their sample were highly 
variable (defined as AS > 50%) and hence derived a 
rate of p = 18 deg ~^. Also in s e arch of highly vari- 
able radio sources, iBecker et al.l (|201GD surveyed the 
Galactic plane and found p = 1.6 deg~^ above a sen- 
sitivity 1 mJy at 4.8 GHz. At low Galactic latitude 
and sampling a m ultitude of timescales (days to years) 
lOfek et akl (|201lh used the VLA to study transient and 
variables in a field covering 2.66 deg^ to a sensitivity 
limit of ~ lOOpJy at 5 GHz. One transient candidate 
was reported and a thorough discussion of the variable 
rad io sources within the fiel d was given. 

iMatsumura et all (|2009l l summarise the detections 
of nine candidate tra nsient sources from th e Nasu 1.4 GHz 
wide-field survey (seelKunivoshi et al.|200 7: Matsu mura et 
20071: iNiinuma et al.l2007l:lKida et al.l2008i:,Nfinuma et all 



interpreted. For Figure [3] we have adopted the conser- 
vative snapshot rate calculated bv lFrail et al.l ||2012|) . 
which assumes no detections. See iFrail et all (|2012l ') 

for further discussion. 

The Allen Telescope Array (ATA: IWelch et al.l2009l ') 
is the first of the next generation telescopes to have 
been used for large area transient surveys. The Pi 
GHz Sky Survey (PiGSS) survey was carried out on 
the ATA at 3.1 GHz and covered 10 000 deg^ of the 
sky with a minimum of two ep ochs for each field w ith 
an RMS sensitivity of ~1 mJy (jBower et al.ll2010l '). 

The first data release (PiGSS-I) is a 10 deg'^ region 
in the Bootes constella tion, with 75 daily observations 
over a 4 month period. iBower et al.l (|2010l ') identify one 
object present in PiGSS that is not in archival radio 
catalogues (NVSS and WENSS). From this they set an 
upper limit on the snapshot rate of transients of p < 1 
deg~^ at 1 mJy sensiti vity at 3.1 GHz. The second 
data release (PiGSS-II; H ower et al]|201lh focused on 
studying the daily and monthly transient and variable 
nature of the field. The upper limits of p < 0.025 
deg~^ at 15 mJy for one day timescales and p < 0.18 
deg~^ at 5 mJy for one month timescales were placed 
on the rates of transient events at 3.1 GHz. 

The other major ATA survey is the The Allen Tele- 

scope Array Twenty-Centimeter Survey (ATATS: ICroft et al.l 
20101') ■ a 690 deg survey with 12 epochs of each field. 



Croft et all (|2010D presented results from comparing a 



20091 . for full details). The sources were detected in a 
drift scan mode, and have flux densities greater than 
1 Jy with typical timescales of minutes to days. (Note, 
however, the questions raised about these results by 
ICroft et al ] |2010l . I2OIII and the discussion about the 

resulting transient event rates bv lOfek et al.|[201ll.) 

As briefly discussed earlier in Section lXTliBower et all 



(|2007l l reported the detection of ten transient radio 
sources in archiv al VLA data at 4.8 and 8.4 GHz. 
iFrail et aP (|2012l ') recently reanalysed this dataset and 
reported that more than half of these transients were 
either caused by rare data reduction artifacts, or that 
the detections had a lower S/N after re-reduction. The 
iBower et al.l (|2007l ') study initially predicted a rate of 
events over the whole sky of p ~ 1.5 deg^^, which is 
diminished to potentially an order of magnitude less, 
depending on how the lower S/N level transients are 



single combined ATATS image with the NVSS. They 
did not detect any transients, placing an upper limit 
of p < 0.004 deg^^ on the transient rate for sources 
with fl ux greater than 40 mJy at 1.4 GHz. ICroft et al.l 
(|201ll ') extended this analysis presenting results for an 
individual epoch-to-epoch comparison for the ATATS 
dataset. They derived a limit of p < 6 x 10~* deg^^ 
for transients brighter than 350 mJy at 1.4 GHz. 

3.4 Low Frequency Surveys 

Since many significant discoveries of radio transients 
have been made at low frequencies, the low frequency 
SKA pathfinders have enormous discovery potential. 

The largest blind survey f or slow transients a t low 
frequencies was conducted bv lLazio et"al] (|2010f ) with 
aOie Long Wavelength Demonstrator Array (LWDA). 
The LWDA is a 16 dipole phased array with all-sky 
imaging capabilities, operating at 73.8 MHz. A total 
of 106 hours of data (with a typical integration time 
of 5 minutes) was searched for transients. No detec- 
tions were made a bove a flux density limit of 500 Jy. 
iLazio et al.l (|2010l ) quote an upper limit on the rate 
of 1 X 10~^ yr~^deg~^ which converts to a rate of 
p < 9.5 X 10~* deg~^ assuming 5 minute integrations. 

In comparison, one of the deepest l ow frequency 
studies to date has been conducted by jJaeger et al.l 
(|2012l ). They used six VLA observations centred on 
the SWIRE Lockman Field at 325 MHz to study sources 
brighter than 0.2 mJy within a field of view of 6.5 
deg^. The timescales of observation ranged from 1 
day to three months. One transient was reported and 
the authors argue that it is most likely coherent emis- 
sion from a stellar flare. Using this one detection, 
jJaeger et al.l (|2012l ') set an upper limit on the snap- 
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shot rate of p < 0.12 deg ^ for sources brighter than 
2.1 mjy at 325 MHz. 

Between the two comprehen sive low frequency tran- 
sient surveys disc ussed above (jLazio et al.l I2OIOI and 
IJaeeer et al] 120121 ') lies a large amount of parameter 
space both in frequency and time to explore. This will 

be made possible by two ne w instruments : the Mu rchi- 

son W idefield Array (MWA: lLonsdale et al . 2 009,: Tingav et 
2012 ) and the Low Frequency Array ("LOFAR: lHeald et al.1 
201l|), both currently being constructed and commis- 



sioned. 



4 Survey Strategy 

The detection and characterisation of radio transients 
requires both the sampling of new regions of param- 
eter space, and the ability to monitor and follow up 
known or newly detected sources. As such, a success- 
ful transient instrument requires a large field-of-view, 
rapid survey capabilities, high sensitivity and dynamic 
range, as well as sufficient angular resolution for unam- 
biguous source identification and follow-up. ASKAP 
will fulfill these criteria with its ability to achieve sub- 
millijansky sensitivity across the entire visible sky in a 
single day of observing. 

Sources that vary on rapid timescales are necessar- 
ily compact. The 10" resolution of ASKAP will permit 
source localisation to ~arcsecond positional accuracy, 
enabling follow-up observations at other wave-bands. 

The key to the efficient detection and monitoring of 
transient and variable sources is repeated observation 
of each field to sample a full range of cadences corre- 
sponding to the natural timescales of variability in the 
underlying source population. In addition to observ- 
ing commensally with other ASKAP programs, we will 
run a campaign of complementary surveys designed to 
achieve our science goals. The VAST surveys build 
on the standard ASKAP imaging pipeline and thus do 
not include other possible m odes of operation suc h as 
"Fly's Eye" monitoring (e.g.. ISiemion et al]|2012l '). 

4.1 VAST Surveys 

In this section we outline our strategy for each of the 
VAST surveys given in Table [2] We have based this 
survey strategy on an estimate of ~ 8500 hours of ob- 
serving time over the first five years of ASKAP opera- 
tions, in addition to continuous commensal observing 
with other ASKAP projects. The survey parameters 
are summarised in Table [J] All of the dedicated (non- 
commensal) VAST surveys will be conducted with an 
observing frequency range of 1130-1430 MHz. 

4.1.1 VAST Wide 

To detect rare events such as GRBs and SNe, and to 
map out the gas density distribution of the ISM across 
the whole sky, a large survey area is required. To 
fully sample the ESE light curves, a regular (prefer- 
ably daily) cadence is desirable. Two years of data are 
required to discriminate between an annual cycle and 
intrinsic variability for IDV sources. 



Hence our planned survey will cover 10 000 square 
degrees (comprised of 400 pointings, each of 30 square 
degrees). Integrating for 40 seconds per pointing will 
give an expected RMS sensitivity of 0.5 mJy beam~^, 
with a total observing time (including overheads) of 6 
hours per day. This survey will be repeated daily for 2 
years, with a total time of 4380 hours. Figure [3] shows 
^kt it will be a substantial improvement over existing 
blind gigahertz surveys. 

4.1.2 VAST Deep 

In our Deep survey we plan to observe 10 000 square 
degrees, with 400 pointings each observed for 1 hour. 
This would achieve an expected RMS sensitivity of 
50 /iJy beam^^. Each of the fields will be observed 
by the EMU colla boration (see Sect ion 14. 2p for their 
continuum survey orris et al ] |201ll ). While both the 
Deep and the Wide surveys observe the same amount 
of sky, they are complementary in cadence and sensi- 
tivity. 

One field, coinciding with a field observed by the 
EMU collaboration, will be repeated daily for one year 
(Deep Single field). The other fields will be repeated 8 
times at irregularly spaced intervals (Deep Multi-field). 
This is the optim al survey to detec t as yet unknown 
source classes fsee lBower et al.ll2007l ) and detect GRBs 
and SNe to larger distances. 

4.1.3 VAST Galactic 

A third component of VAST is a survey of 750 square 
degrees of the Galactic plane. This is a total of 30 
pointings, plus an additional 3 pointings to cover the 
Large and Small Magellanic Clouds. Each pointing 
will be observed for 16 minutes, achieving an expected 
RMS sensitivity of 0.1 mJy beam~^. The Galactic 
plane survey will require 9 hours of observing in a sin- 
gle day, repeated weekly for 1 year and supplemented 
by additional epochs at longer intervals. 

4.2 Commensal Observing 

In addition to the planned surveys, we will operate our 
transient detection pipeline commensally on all other 
ASKAP observations. Two pr ojects, EMU (The Evo- 
lutionary Map of the Universe : I Norris et al.ll201ll fland 
WALLABY (The ASKAP HI All-Sky Survey; Korib- 
alski, et al., in preparation)0, are particularly relevant 
as they will conduct extensive (possibly commensal) 
wide-field surveys: 

EMU: a deep (10 ^Jy beam"^ RMS) survey of the 
entire Southern Hemisphere, to a declination of 
5 = -f30°. 

WALLABY: a survey of neutral Hydrogen covering 
two thirds of the sky to a redshift of z — 0.5. 

These projects plan to operate in 'point and shoot' 
mode, covering some region of the sky with signifi- 
cant integration time (~ 12 hours) per pointing and 
no revisits. This provides an excellent opportunity for 



^Scc http://www.askap.org/cmu 
'^See ^http: / /www.askap.org/waIlaby[ 
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detection and monitoring of transient and variable be- 
haviour on timescales of less than a day. 



4.3 VAST Predictions 

In Figure [3] we show the predictions parameterised via 
the method described in Section [S] using the VAST 
survey parameters given in Table[3] Note, the horizon- 
tal lines show the upper limits that would be placed 
if no transient or variable sources were detected. The 
vertical lines show the RMS that would be expected 
for the final integration time per image (see RMS sen- 
sitivity in Table[3|. For the commensal survey the best 
RMS is found after six hours (not 12), assuming that 
two back-to-back six hour images are needed to search 
for transients, i.e., that there are no repeat observa- 
tions. It can be seen from Figure |3] that the VAST 
surveys make orders of magnitude improvements on 
previous surveys. 

VAST- Wide will use 4320 hours of dedicated ob- 
serving time, which is the greatest of all the VAST 
sub-surveys (excluding commensal observing). It will 
observe 10 000 deg'^ for 40 seconds every day for two 
years, possibly longer. Observing for 40 seconds with 
ASKAP yields an RMS of ~0.5 mjy beam"^ and pro- 
duces eight five-second images per observation. This 
is represented on Figure [3] and it casts VAST- Wide as 
less competitive than, for example, VAST-Deep (both 
multi-field and single field). The limits shown in Fig- 
ure 13] for VAST- Wide however do not include upper 
limits on the rate and RMS calculations via combin- 
ing data from multiple days observing. 

For example, all the 40 second images taken daily 
for an entire week could be combined together and 
compared and searched for transients with the com- 
bined images from the next week, and so on. If tech- 
nically feasible, such an approach will significantly in- 
crease the competitiveness of the VAST- Wide survey. 
Figure [3] also does not convey how many known radio 
sources will be monitored in each case. VAST- Wide 
will monitor all known sources in the Southern hemi- 
sphere above 2.5 mJy (5a") every day for two years, 
whereas VAST-Deep (one field) will monitor only a 
subset (30 deg^) to a greater sensitivity of 250/iJy 
{5a). We can approximate the ratio of sources that 
wiU be observed with VAST-Deep {Noeep) to VAST- 
Wide (Nwide) using the following expression (assum- 
ing a Euclidean source population): 



D eep 



A 



Deep 



Nwii 



A 



Wide 



a Deep 
O Wide 



10% (6) 



5 Data Processing 

The success of VAST hinges on our ability to detect 
transient sources on the timescales on which images are 
produced by ASKAP, at the same time as monitoring 
known variable sources that are in our field. Firstly, 
how to effectively deal with the large volumes of data, 
and secondly, how to find the 'needles in a haystack' 
— the most interesting transient sources. As discussed 
in Section [31 current blind surveys at 1.4 GHz sug- 
gest that less than 0.5% of sources ob served by VAST 
will b e transient or highly variable (|Bannister et al.l 
l2011al ). 

The operation of the ASKAP data pipeline is de^ 

scribed in the Science Processing Document ( Cornw cU et al.l 
|2012i). The following sections address some of the 
issues most relevant to VAST, and discuss the de- 
sign and operation of the VAST transient detection 
pipeline. The overall functionality of the VAST pipeline 
is illustrated in Figure [4] and described in Table (5) A 
discussio n of the p ipeline design is also provided by 
[Banvcr (|2011| '). 

5.1 Data Rates 

ASKAP is designed to be a real-time telescope, with 
most of the data processing done in online mode. Data 
from the phased array feed elements in each antenna 
will be sent to the beam former to create primary 
beams. The primary beams from each antenna are sent 
to the correlator to create visibility data. These data 
are then sent to the central processor for processing 
into science-ready data products: images and spectral 
cubes. The data rate from the correlator to the cen- 
tral processor is expected to be 2.5 GB/s, generating 
about 200 TB of raw data every 24 hours. 

The central processor will operate three data pro- 
cessing pipelines in parallel: one pipeline will produce 
data cubes designed for spectral line analysis, with 
high frequency resolution; a second will produce data 
cubes designed for continuum surveys, with full po- 
larisation and higher angular resolution, but reduced 
spectral resolution; the third pipeline is the transient 
detection pipeline which will produce data cubes with 
lower frequency resolution, but operating on 5-second 
timescales (see Table [4]). Each 5-second image cube 
will be up to 16 GB (depending on the number of fre- 
quency channels generated and stored), which means 
that the transients pipeline will process around 12 TB 
of data an hour. Retaining some coarse frequency 
resolution will enable post-correlation rejection of ra- 
dio frequency interference (see section [53] below) and 
also provide some minimal spectral information for de- 
tected transient sources. 



where in this example Aoeep is the survey area of 
VAST-Deep (single field, 30 deg^); Awide the area of 
VAST- Wide (10,000 deg^). aoeep and awide represent 
the survey sensitivities respectively. VAST- Wide will 
therefore monitor large numbers of bright transients 
and variables, such as IDVs and AGN. VAST-Deep 
will be better optimised for studying the faint end of 
the transient and variable population. 



5.2 Transient Detection Approach 

Source detection is central to all upcoming radio tran- 
sient surveys, and also those in other wave-bands such 
as the Large Synoptic Survey Telescope transients sur- 
vey. In the image domain, transient source detection 
has typically followed one of two approaches: image 
subtraction or catalogue-based detection. 
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Tabic 4: VAST Image Products 



Polarisations 


1-4 


Frequency channels 


30 


Field of view 


7.5 deg2 


Angular resolution 


10" 


Integration time 


5 seconds 


Image size (pixels) 


4096 X 4096 


Image cube size 


up to 16 GB 
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Supernova search projects and other optical as- 
tronomy projects have typicall y used the image sub- 

tract ion approach (for example. I Riess et al.ll2004 : lAlcock et 
120001 ). In its simplest form, this involves subtracting 
the current image from some higher sensitivity mas- 
ter image. In practice, some corrections are required, 
for example each image can be convolved with a ker- 
nel to account for coordinate shif ts and distortions 
(|Alard fc Luptonlll998l : lAlardll2000D . 

In radio images, the incomplete u-v coverage means 
that image subtraction often results in artefacts that 
can then be mistakenly identified as transient sources. 
A more robust approach is to firstly extract and mea- 
sure sources in each image, and then search for tran- 
sients by catalogue cross-matching and lightcurve anal- 
ysis. The approach taken in the VAST transient de- 
tection pipeline is base d on techniques described by 
iBannister et al.1 (|2011al l. developed for analysing the 
Molonglo Observatory Synthesis Telescope archive. This 
is also similar to the approaches t aken in the Allen 

Teles cope Array transients projects (jCroft et aI.ll2O10l. 

20111 ) and the LOFAR transient detection pipeline |Swinbani3 
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20071 ) 



We are also investigating the techni que of tran- 
sient detection in the vi sibility domain (jTrott et al.l 
I2OIII : iLaw fc Boweill2"012l ) . These approaches have the 
advantage of dealing with well-behaved noise prop- 
erties of visibilities and avoiding artefact s caused by 
the im aging and deconvolution processes. iTrott et al.l 
(I2OI1I ) use statistical decision theory to derive detec- 
tion hmits of visibilities using models of the source and 
calibr ation effects of interferometers. iLaw fc Bowed 
(|2012l ) present a computationally-simple, calibration- 
independent algorithm for transient detection based 
on the bispectrum, an interferometric closure quan- 
tity. By subtracting visibilities in time, the bispectrum 
can be used to quantify the significance of transients 
throughout the field of view. Visibility-based algo- 
rithms have great potential in the study of transients, 
because they take advantage of a priori information to 
reduce computational burden. This makes real-time 
processing easier, which may be a critical limitation 
of the large, computationally-demanding surveys de- 
scribed here. More work is needed to understand how 
these algorithms work on a large scale and whether 
they can be implemented for real-time analysis. 



gure 4: Illustration of the VAST pipeline func- 
tionality. A description of each stage is given in 
Table [3 



5.3 Rejection of Radio Frequency In- 
terference 

One of the most significant challenges for all radio 
telescopes is the impact of radio frequency interfer- 
ence (RFI) from man-made sources. Although the 
ASKAP site was chosen for its radio-quiet properties, 
signals from satellites used for communications, nav- 
igation, remote sensing and the military pervade all 
areas of the Earth. Similarly aircraft communication 
signals are present even in very remote areas. In addi- 
tion, self-generated RFI is a challenge as faster digital 
equipment is used in the signal pathways. There is 
no universal method of RFI mitigation that is satis- 
factory for all types of observing techniques and RF I 
signals (jFridman fc Baanll200ll : iBriggs fc Kocj|2005l ). 
For the data rates expected with VAST, an automated 
approach with a combination of techniques will be es- 
sential. 

RFI mitigation techniques range in sophistication 
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Table 5: Functionality of the VAST pipeline components, as illustrated in Figure HI 



Function 



Description 



Quality control & RFI 
mitigation 



Source detection 



Source association 



Source monitoring 



Quality control & gain 
calibration 



Light curve creation 



Transient detection 



Multi-wavelength cross- 
matching 



Transient classification 



Quality control & final 
data products 



Initial quality control on the image cubes output from the 
ASKAP pipeline. Poor quality images will be flagged or 
ignored. Possible RFI mitigation strategies are discussed 
in Section [Ql 

Blind source finding is performed on every image as dis- 
cussed in Section 15.41 This produces a list of detected 
sources with position, flux, size and local RMS. 

Each detected source is associated with one from the mas- 
ter list by finding the master source at the same position as 
the detection. If no master source is found at that position 
then a new master source is added. 

Master sources known to be in the field that were not de- 
tected by the source finder are measured by attempting to 
fit a Gaussian at the position where the source is expected 
to be. If no source can be fit then an upper limit on the 
flux at that position is taken instead. 

The image is assessed for quality by comparing the mea- 
sured parameters of the detected sources to a list of well- 
characterised sources u sing the techniques described by 
iBannister et all (|2011al) . 

All measurements for each source are collated into a ra- 
dio light curve — a collection of flux measurements for all 
epochs where the source was observed. 

The updated light curves are analysed to detect transient 
or variable behaviour. These sources will go on list for 
monitoring, and cross-matching with external catalogues. 

Sources of interest are cross-matched with external cat- 
alogues to provide additional information on the likely 
causes of variability, and to collect extra information for 
classification and identification. 

Sources that exhibit transient or variable behaviour are 
classified using a machine-learning algorithm. Possible ap- 
proaches are discussed in Section [5751 

After a final quality control stage, sources that appear gen- 
uine with a high degree of confidence will be released to 
the community via VOEvent notifications as discussed in 
Section [5T6l These alerts, along with science quality cata- 
logues of transient and variable sources will be made public 
through the VAST archive. 



from simple blanking of corrupted data to adaptive 
cancellation based on coherent detection. Simple blank- 
ing works well for strong short duration RFI, but is 
completely ineffective in removing low-level time- variable 
RFI that could mimic transient sources. Blanking re- 
moves signal as well as interference and reduces the 
sensitivity achieved in a given time. There is a growing 
literature on more sophisticated ideas for RFI elimina- 
tion. Some of the techniques that we are exploring for 
VAST include: 

1. Post-detection flagging based on recognition of 



data that are affected by interference. This has 
been the traditional approach in radio interfer- 
ometry and is appropriate for moderately strong 
interference, but it does come with a loss of sen- 
sitivity. 

2. Adaptive c ancellation technique s such as those 
descri b ed inlMitchell et al.1 (12005^ ; iBriggs fc Kocd 
(120051 ): iKesteven et al.1 (|2010l ). These require a 
separate reference antenna or beam on the sky 
to sample the RFI. An algorithm is used that 
adapts filter characteristics to match the RFI 
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spectral footprint. This approach requires a higli 
sensitivity copy of the interferer. 

Spatial filtering techniques that use the loca- 
tion of the known interferer or the relative ar- 
rival times of the signals to identify and sepa- 



nal fe.E. iKocz et all l2010l: 


Leshem et all I2OOOI: 


Ellineson & HamDsonl2002l: 


Smolders & Hamnsonl 


2002t). These are aoDrooriate for stationary or 



slowly moving RFI sources and require less sensi- 
tivity to the interfering source than the adaptive 
cancellation techniques. 

4. Adaptive nulling techniques, which employ the 
capacity of the array to place response nulls on 
the RFI source. These are appropriate for very 
strong interfering sources and for telescopes such 
as ASKAP with independently steerable beams. 

5. Post-correlation filtering techniques use software 
algorithms to remove RFI from interferometric 
data b y subtraction of a weak, persistent RFI 
signal (|Athrevall2009l ^ or by implementation of 
low-pass time and frequency filte rs to remove 
variable RFI (|Qffringa et al.H2012l '). The advan- 
tage of these techniques is that they generally 
have a lower sensitivity requirement for detect- 
ing the interference. However, radio interferom- 
eters may not be able to use this approach to 
mitigate some types of slowly varying RFI sig- 
nals. 

5.4 Source Detection 

Finding and characterising sources in each of the ob- 
served images are critical first steps for transient de- 
tection in the image domain. False positives in source 
finding have the ability to confound the detection of 
transients, reducing the confidence with which a tran- 
sient can be detected. This has been identified as a 
problem in previous surveys, with source-finding errors 
con tributing to the false detection rate (see, for exani- 
ple.lBannister et al.ll2011al : ICroft et al.l201ll : iFrail et all 

l2om 

The source finding packages that are cu rrently avail- 
able (for example the Miriad packa ge sfind ([Hopkins et alJ 
l2002l 'l and the AIPS package VSAD (jCondon et al.llQQST l'l 
are all able to achieve a high degree of completeness 
and reliability, typically better than 95% on standard 
radio continuum images. These packages have been ex- 
tensively used in large-scale radio surveys. However, 
even a small number of missed real sources, or false 
positives, can cause problems in the identification of 
variable and transient sources. 

In recent surveys the best source finding results 
have been obtained by using the above mentioned source 
finding packages in conjunc tion with a pipeline of other 
programs (see, for exam ple, iMauch et al .i2003l : lMurphv et 




l2007l : lNorris et al.ll200^ '). The refinement of the source 
finding pipeline is often done after all observations 
have been completed and the idiosyncrasies of the data 
are well understood. A lot of time and effort is spent 
normalising the input images, and filtering the result- 
ing catalogues. In many cases a subset of the final 



Figure 5: The false detection rate (FDR) 
for each of the sou rce finding algorithms (from 
iHancock et al ] l2012l ). No falsely detected sources 
are expected above an SNR of 5 for the area of sky 
simulated. The non-zero rate of false detections is 
due to poor source characterisation. 



source catalogue is validated 'by eye' to ensure that 
the brighter sources are correctly represented. This 
approach produces high quality catalogues and is well 
suited to single epoch surveys, as is evidenced by the 
high level of completeness and reliability that these 
surveys are able to obtain. 

For VAST, the volume of data will require a fully 
automated source finding process. The inability to 
archive every image that is observed for VAST will 
require that the source finding is done in real-time, 
and without the possibility of reprocessing. The large 
number of sources that are expected for a single day's 
observation means that even with completeness and 
reliability of 99%, a source finding algorithm will be 
finding up to ~ 10, 000 false detections and missing up 
to ~ 10, 000 real sources each day. These false detec- 
tions and missed sources far outweigh the number of 
true transient and variable sources that we expect to 
detect. 

As part of the VAST design study process. iHancock et al.l 
(|2012l l have undertaken a detailed analysis of some of 
the most widely used source finding packages i n astron- 
omy, including the newly developed Sela vy (|Whitinel 
l2012l 'l and Aegean (|Hancock et al.1 120121 ') algorithms. 
For compact sources, we find that the task of source 
finding (distinguishing between source and background 
pixels) is a solved problem that the various programs 
all achieve equally well. However, the task of source 
characterisation (measuring the various parameters of 
ailhe found sources) is only accomplished well in the 
case of single isolated sources. For diffuse sources, the 
task of source finding is equally well achieved, how- 
ever t he task of source ch aracterisation is not well de- 
fined. iHales et al.l (|2012l ) address some of the issues 
of extended source measurement. However, since we 
expect transient sources to be compact, this is not a 
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Figur e 6: The complete ness of each catalogue 
(from [Hancock et al.ll2012l ) as compared to the in- 
put catalogue. The coloured curves represent the 
completeness of the named source finders. The 
black dotted curve represents the expected com- 
pleteness of an ideal source finder. The plateau 
in completeness above and SNR of 10 is due to 
sources that are poorly characterised. 



significant issue for VAST. 

Figures[5]and[6]show the false detection rate (FDR) 
and completeness of each of the source finding algo- 
rithms. The FDR and completeness were measured 
by injecting a known population of sources into a sim- 
ulated image and processing this image with each of 
the different source finding algorithms. A comparison 
between the input and output catalogues resulted in: 
a set of sources common to both (real detections), a 
set of sources in the input but not the output (non- 
detections), and a set of sources that were in the out- 
put but not the input (false detections). The set of 
falsely detected sources consisted entirely of sources 
within an island of pixels that contains multiple com- 
ponents. These sources were so poorly characterised 
that they were not able to be linked to a source within 
the input catalogue. Sources that are poorly char- 
acterised also reduce the catalogue completeness. In 
Figures [S] and El the main reason that the source find- 
ing algorithms have non-zero FDR and less than ideal 
completeness is due to poor source characterisation, 
rather than poor source finding. In order to better 
characterise islands of pixels that contain more than 
a single component the Aegean program analyses the 
curvature of an image. As can be seen in Figures [5] 
and[6l Aegean is able to characterise islands of pixels 
that contain multiple components better than any of 
the other source finding approaches, typically achiev- 
ing a completeness of over 99.9%. Further details of the 
source finding analysis and the operation of the Aegean 
source finding algorithm are presented bv lHancock et akl 
(|2012l ). 



5.5 Classification of Transients 

Given the large number of sources that we will likely 
detect with VAST, prompt classification will be critical 
for early identification and multi-wavelength follow- 
up. This is a challenging problem due to the sparsity 
and heterogeneity of the available data. There are two 
sets of information that can be used for this task: the 
light curve itself, and the existing archival information 
about that source or location on the sky. 

Light-curve classification has been successfully used 
to classify variable stars at optical wavelengths (e.g. 
[Richards et"al]|201ll . l2012l l . In order to explore the fea- 
sibility of light-curve classification for VAST, we used 
several approaches. 

Firstly, we simulated light curves of five transient 
source types: SNe, RS CVn flare stars (fStar_RSCVn), 
M-class flare stars (fStar_dMe), XRBs, novae; two vari- 
able source types: IDVs and ESEs and a non-varying 
background (BG) for comparison (Lo et al. in prepara- 
tion). We experimented with light-curve classification 
on the simulated light curves in two scenarios: 

• archival classification where the classification is 
performed after a significant amount of light curve 
data has been collected; and 

• online classification, which is classification of par- 
tial light curves in real-time for the purpose of 
triggering appropriate follow-up observations. 

Secondly, we used light curves from the Palomar-Quest 
and Catalina Real-Time Transient Surveys to experi- 
ment further with a range of machine learning meth- 
ods. The results of these experiments are discussed in 
the sections below. 

5.5.1 Archival Classification 

We approached the classification problem using the 
standard methodology of supervised machine learn- 
ing with 10-fold cross validation. First we extracted 
a set of features from each light-curve. These include: 
time domain features which are simply the flux mea- 
surements at each time point (tme), statistical features 
(stat) suc h as the mean, s tanda rd deviation and skew, 
following [Richards et al.l (|201lD : statistical (stat) fea- 
tures built cumulatively {stat- cum) after observing a 
fixed number of observations in multiples of k (k, 2k, 
etc.); periodic features derived from the Lomb-Scargle 
periodogram {Isp); and wavelet features derived from 
the Haar wavelet (wlet). We then created two more 
feature sets that are concatenations of the above. Fea- 
ture set all-reps is the concatenation of stat-cum, Isp, 
and wlet, and feature set all concatenates tme, stat- 
cum, Isp, an d wlet. Details of all s e ven fe ature sets are 
discussed bv lRebbapragada et al.l (|2012l '). 

We used these features as input to a set of ma- 
chine learning algorithms, which, when given a train- 
ing set of examples, will learn the discriminative fea- 
tures for each so urce type. Using the Weka toolkit 
(jHall et al.ll2009l ). we experimented with three classi- 
fication algorithms: S upport Vector Ma chines (SVM; 
ICortes fc VapnikI 19951 1. Random Forests ('Brciman"20"oil) 
and the J48 Decision Tree (fQuinlan. .IDSS l . Figure [7] 
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Figure 7: Comparison of overall classification pt 
formance (accuracy percentage) by feature set 
{tme, stat, stat-cum, wlet, Isp, all-reps, and alt) 
and classification algorithm (J48, Random Forest, 
and SVM), with results grouped by classification 
algorithm. 



shows the classification accuracy in the archival set- 
ting. The Random Forest is the best performing classi- 
fier when used with the combined feature sets, achiev- 
ing 93% accuracy. This suggests archival classification 
using light curves alone is feasible with a reasonable 
degree of accuracy, given the large number of sources. 
However, the performance on real light curves is likely 
to be lower since our simulated light curves have ideal 
Gaussian noise and are aligned to the starting time of 
the event. 

5.5.2 Online Classification 

The aim of online classification is to identify sources 
of interest as soon as possible after new observations 
arrive. The earlier a classification is made, the more 
likely that a follow-up observation can be scheduled, 
and so online classification should run continuously as 
part of the real-time transient detection pipeline. We 
evaluated online classification using the first 30 days 
of simulated light curves. With so few observations, it 
is not possible to extract periodogram or wavelet fea- 
tures, hence we were limited to the statistical feature 
set and the time series values themselves. 

Figure |8] shows the confusion matrix for our best 
performing algorithm: the J48 classifier built with stat 
features extracted from 30 days of observations. The 
overall classification accuracy achieved was 50% across 
eight source types. This is better than the baseline 
expected from random classification (12% accuracy) 
but probably not sufficient to be useful in schedul- 
ing follow-up observations. Accuracy was improved by 
grouping source types with similar light-curve charac- 
teristics together (grouping SNe and Novae, BG and 
ESE yielded an accuracy of 70%). Finally, accumulat- 
ing more data led to higher accuracy {r^ 80% accuracy 
with 100 days of observations). 
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Figure 8: Online classification confusion matrix 
shown as a heat map (red = higher accuracy, blue 
= lower accuracy). The x-axis is the classified out- 
put class and the y-axis is the actual input class. 
We used 200 sources of each type and the num- 
ber in each cell represents how these sources are 
classified. Overall accuracy at 30 days is 50%. 



5.5.3 Other Machine Learning Approaches 

Considerable insight into the classification problem can 
be gained from the current efforts on cl assification of 
transients in optical synoptic sky surveys (iBloom et al.l 



., survevsj i 

200s ;|Don^ek et al.ll2008l:[Mahabal et al.l2008J bl. l2010l . 



201ll : lDiorgovski et al.ll2011^ . We experimented with 
several approaches, using the data from the Palomar- 
Quest and Catalina Real-Time Transient Surveys, as 
well as selected data sets from the literature. 

Similar to the work described in Section 15.5.11 we 
parameterised observed optical light curves as feature 
vectors, and applied various machine learning tech- 
niques (e.g., decision trees) to their classification. A set 
of decision trees were trained using the feature vectors 
for various combinations of classes. To reduce the di- 
mensionality of the input space, we applied a forward 
feature selection strategy that consisted of selecting 
a subset of features from the training set that best 
predict the test data by sequentially selecting features 
until there is no improvement in prediction. Each tree 
was built using the Gini diversity index as criterion for 
choosing the split; the splitting stopped when there is 
no further gain that can be made. To avoid over-fitting 
we used a 10-fold cross validation approach. The deci- 
sion trees were pruned in order to choose the simplest 
one within one standard error of the minimum. When 
tested on a data set consisting of light curves of blazars, 
cataclysmic variables, and RR Lyrae stars, we achieved 
completeness in the range of ~ 83% - 97%, and con- 
tamination in the range of ~ 4% to 13% (Donalek et 
al., in prep.). In future work we will apply this to the 
classification of radio light curves. 

Another approach we experimented with used two- 
dimensional distributions of magnitude changes for dif- 



www.pu blish . csiro. au /journals/pasa 



23 



ferent tim e baselines for all possib le epoch pairs in the 
data set (|Diorgovski et ai]|201ll . Moghaddam et al., 
in preparation). These two-dimensional (Am, At) his- 
tograms can be viewed as probabilistic structure func- 
tions for the light curves of different types. Template 
distributions for different kinds of transients and vari- 
ables are constructed using the reliably classified data 
with the same survey cadences, S/N, etc. For any 
newly detected variable or a transient, corresponding 
(Am, At) histograms are accumulated as the new data 
arrive, and we used a variety of metrics to compute 
the effective probabilistic distances from different tem- 
plates. Initial testing indicates that we can obtain a 
high classification accuracy using this approach. We 
plan to generalise it to include triplets or even higher 
order sets of data points for multi-dimensional his- 
tograms. 

One important lesson so far is that existing archival 
and contextual data will play a critical role in classi- 
fication. This includes a spatial context (i.e., what 
is near the observed event on the sky — a possible 
host galaxy, a cluster, a SN remnant, etc.), the multi- 
wavelength context (has it been detected at other wave- 
lengths, what is its broad spectral energy distribu- 
tion), and the temporal context (what was its flux 
variability or a detection history in the previously ob- 
tained data). Some of this information can be readily 
extracted from the archives (e.g., flux measurements 
from different wavelengths, light curves at that loca- 
tion), but some - spatial context in particular - require 
a human judgment, e.g., is the apparent proximity to 
other objects (galaxies, clusters, etc.) likely to be rele- 
vant, and if so, what does it imply about the transient? 
Human inspection of vast numbers of transients does 
not scale to the massive data streams such as those 
contemplated here. We are currently experimenting 
with crowdsourcing approaches to harvesting of the 
relevant human pattern recognition skills and domain 
expertise, and their translation to machine-processable 
algorithms. 

5.6 Alerts and Follow-up 

ASKAP policy specifles that all data and results must 
be released as soon as they have passed scientiflc qual- 
ity control. This is particularly important for tran- 
sient detections, in which multi-wavelength follow-up 
is critical. VAST will generate publicly available data 
products on the following timescales: 

Seconds: Triggers for individual transient events will 
be released through services such as VOEvent- 
NelQ as soon as they are detected. These events 
will have quality flags representing our confl- 
dence about the detection and will allow imme- 
diate follow-up of new discoveries by the wider 
community. 

Minutes: Triggers with value-added data such as clas- 
siflcations based on multi-wavelength archival 
crossmatching will be released as soon as our 
source classiflcation pipeline has run. 

^See |http://voeventnet. caltech.edu [ 



Days: Measurements of every source in every field will 
be entered into our lightcurve database in near 
real-time. This will be available immediately to 
the VAST collaboration. Once quality control 
has been completed, updates to the database 
will be made publicly available. 

Months: Scientific quality catalogues will be constructed 
from our lightcurve database. The results will 
be published and data made public as soon as 
possible after analysis is complete. 

The quality control process for alerts will consist of 
a multi-tiered system in which there is automatic fil- 
tering of a large fraction of false detections followed 
by manual inspection of the remaining sources. This 
approach is taken by several existing transients sur- 
veys, for example in the Palomar Transients Factory 
(|Gal-Yam et al.ll201lf ). The detailed implementation 
of this general approach will be developed and refined 
during commissioning. 

In addition, all images from VAST will be made 
available through the standard ASKAP archive. Where 
possible the VAST pipeline will use standard Virtual 
Observatory protocols such as VOEventfl 

5.6.1 Very Long Baseline Interferometry 

As noted earlier, radio sources that are highly variable 
on short timescales are necessarily intrinsically com- 
pact. ASKAP's angular resolution of ~10"will be high 
enough to localise variable sources on the sky and mea- 
sure variability, but will be orders of magnitude too 
low to resolve the evolving structures that give rise to 
the variations in integrated flux. In order to resolve 
the structure and evolution of the radio source, much 
higher angular resolution imaging is required, provided 
by Very Long Baseline Interferometry (VLBI). 

Alerts from the VAST pipelines will be used to trig- 
ger observations using the Australian Long Baseline 
Array (LB A), now including telescopes from Western 
Australia (ASKAP) to New Zealand (Warkworth), al- 
lowing milliarcsecond-scale imaging. Due to the very 
large number of VAST alerts, targets for VLBI follow- 
up will have to be chosen selectively and will need to 
be restricted to the most interesting or rare classes of 
transient and variable sources. The Australian LBA 
is very well suited to this type of follow-up, as evi- 
denced by previous programs that have taken triggers 
for transient radio s ources from MOST or t he ATCA 
(|Tingav et al.lll995l : [Miller- Jones et aLlbOld l. 

The current LBA has a sensitivity of 60/xJy beam~^ 
for a one hour integration at 2.3 and 8 GHz. This is 
comparable to the sensitivity of ASKAP at 1.4 GHz. 
With ASKAP incorporated as an element of the LBA, 
the sensitivity at 1.6 GHz would be 40/iJy beam~^ for 
a one hour integration. Hence a 12 hour observation 
with the LBA would be able to detect any unresolved 
source found by the VAST surveys. 

The recent advent of real-time eVLBI on the LBA 
means that the VLBI observations can be rapidly pro- 
cessed and analysed. Thus subsequent observations 

^See ptp:// www.ivoa.net/Documents/VOEvent [ 
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can be adapted based on the initial results. This facil- 
ity will benefit studies of new transients that will be 
detected by ASKAP. 



6 Conclusions 

Next generation radio telescopes such as the EVLA, 
LOFAR, MWA and ASKAP will have a substantial im- 
pact on our ability to discover and monitor radio tran- 
sient and variable sources. The VAST survey will play 
an important role in the study of a range of phenom- 
ena, from GRB afterglows and supernovae through to 
cataclysmic variables and radio flare stars, as well as 
propagation effects in the interstellar medium. 

The three VAST surveys (VAST- Wide, VAST-Deep 
and VAST-Galactic) will allow us to do a comprehen- 
sive search of radio transient parameter space. It is 
likely that in addition to discovering a large number 
of known classes of objects, we will also discover new 
types of transient and variable phenomena. Our auto- 
matic transient detection pipeline will be able to run 
simultaneously with other major ASKAP observing 
projects, meaning that we will also analyse data in 
commensal mode. 

All data products produced by VAST will be re- 
leased to the wider community as soon as possible after 
passing quality control standards. Wherever possible 
this will be done through the use of Virtual Observa- 
tory tools and protocols. 

In this paper we have summarised the science goals 
of VAST, described our current results in algorithm de- 
velopment for source-finding and transient light curve 
classification, and made predictions about the likely 
source populations that will be observed by VAST. 
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